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I. INTRODUCTION 


Coupling between the atmospheric circulation and the earth '.s magnetic 
field is strongly suggested by the evidence presented in the literature. 

This evidence covers a wide spectrum of space and time scales, is usually 
in the form of correlations, and has been given for all levels of the 
atmosphere. For example, Flohn (1952) showed that the meteorological 
equator (the ITCZ) is more nearly parallel to the geomagnetic than to the 
geographic equator and that the polar- vortex at 200 mb is more nearly 
centered on the geomagnetic than on the geographic pole. King (1974) has 
shown- that the isolines of total ozone are similar to the- isolines of 
magnetic field strength; and Belmont, at al, (1974b)', showed that the 
contours of the amplitude of the semiannual wave in zonal wind at 50 km 
are more nearly congruent with the geomagnetic, rather than geographic, 
coordinate system. The mechanisms which give rise to these correlations 
are not yet fully understood. 

It should be determined whether the atmosphere or the geomagnetic 
field, or neither, is the independent variable responsible for correlations 
such as the above. If the geomagnetic field is the independent variable 
for a given relationship, then meteorologists ought to include that rela- 
tionship in their studies. For example, a high latitude source of NO 
produced by cosmic rays , which enter the atmosphere at latitudes determined 
by their interaction with the geomagnetic field, is now being included in 
studies of the ozone budget (Crutzen, et al., 1975). On the other hand, if 
the atmosphere is found to be the independent variable for a given relation- 
ship, then such knowledge may be useful to space scientists, but meteorolo- 
gists need not consider that geomagnetic relationship in studies of the 
atmospheric circulation. 

The purpose of the present study is to find relationships between 
stratospheric parameters, 30-60 km, and geomagnetic field parameters. The 
mid- and upper-stratosphere may respond dramatically to geophysical events 
(e.g., Zadvernjruk, 1973), and by studying relationships at high altitudes 
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it may be possible to more readily identify some coupling mechanisms between 
the atmosphere and geomagnetic field. Hopefully, this could help explain 
which are the independent variables for some of the known relationships 
between geomagnetic and meteorological data. The method used here will be 
to compare temporal variations of wind and temperature at rocketsonde 
stations (MRN data) with time variations of the vertical (Z) and horizontal 
(H) geomagnetic field intensity at a nearby geomagnetic observatory. 

The frequency range of time variations which can be studied is limited 
only by the time distribution of MEN data, as the geomagnetic data are 
taken hourly in a (usually) continuous sample. The MRN data are sufficiently 
plentiful to define variations longer than a month, so a major portion of 
the study deals with periodic analysis, the quasi-biennial oscillation, and the 
first three harmonics of the annual wave. On time scales of a season or 
less, midwinter sudden stratospheric warmings are the most spectacular 
events. While the MRN data are too sparse to perform case studies of 
individual warmings, it is possible to stratify all years according to 
whether or not a major warming occurred. This procedure has been used to 
study differences of the variance spectra in MRN and geomagnetic data during 
years when major warmings occurred compared with the other years. Finally, 
results are given for a superposed epoch study of the changes in strato- 
spheric temperature over a few hours time at Fort Churchill using solar 
sector boundary crossing dates as the key events . 

II. DATA 


Meteorological rocket (MEN) data 1960-72 were obtained from the World 
Data Center, Asheville. Station locations and the nearest geomagnetic 
observatories are given in Table 1. Further details concerning the MRN 
data and results of periodic analysis of wind and temperature have been 
given in Belmont, et al. (1974a), and Nastrom and Belmont (1975). 

Daily mean values of the geomagnetic field elements for years corres- 
ponding to the MRN data were obtained from the World Data Center, Boulder. 
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Observatories used in this study are listed in Table 2. Generally, geomag- 
netic data prior to 1960 vere not used in order to make the periods of 
record of the MRN and geomagnetic data as compatible as possible. Also, 
the analysis was limited to observatories near a MRN station (see Table 1). 

At most observatories the field elements given are declination (D) , 
horizontal (H) , and vertical (Z) field intensity. At the Canadian stations 
indicated in Table 2 they were reported as X, Y, and Z; but daily values' 
were converted to D, H, and Z prior to further processing. The observatory 
at San Juan was moved at the end of 1964, and it was necessary to adjust 
the base .'line of 1960-1964 data to be consistent with 1965-1972 data. Also, 

i 

the observatory at Honolulu was moved during 1960, so data for I960 were 
not used there. 

III. RESULTS 

A, PERIODIC ANALYSIS OF GEOMAGNETIC FIELD ELEMENTS 

1. Procedure 

Significant peaks at 12 and 6 months are present in the H and Z spectra 
(Currie, 1966), and resolving them as discrete lines provides the possibility 
of studying their phases as well as amplitudes. However, secular trends are 
often very large in the H and Z data, and failure to remove them prior to 
periodic analysis can lead to inconclusive results (Chapman and Bartels, 

1940, Chapter 16; Currie, 1966). Inspection of plots of our time series for 
1960-1972 (not shown here; see Chapman and Bartels, 1940, p. 132) indicated 
that a parabola can be used to effectively remove the secular trend. This 
technique -is more desirable than other filters because no data are lost at 
each end of the time series. Although a parabolic trend line does interact 
with the approximately 11 year cycle found in H and Z, numerical tests 
made using synthetic time series show that the error in amplitude and phase 
of the 11 year cycle, and shorter periods, is less than 4% after removing a 
parabolic trend from a time series 13 years long. 


Time series of mean daily H and Z values are characterized by a 
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relatively steady background which may be dramatically interrupted during 
a geomagnetic disturbance. The effect of geomagnetic storms, which last 
for only hours or days, could be thought of as a very large amplitude 
high-frequency variation which occurs more during some years than during 
others . As the purpose here is to study the month-to-month changes of 
the background geomagnetic field, it is desirable to remove the aliasing 
of monthly data caused by the irregular occurrence of geomagnetic storms. 
One method to achieve this is to use only non- disturbed days when comput- 
ing monthly means . This method has the drawback that disturbed days can 
be identified only on a subjective basis. If the disturbed days form 
only a small part of total daily values in a month, however, then an 
objective and nearly as effective method is to use monthly medians rather 
than monthly means. Periodic analyses were made using both monthly 
median and monthly mean data. The resulting amplitudes and phases differed 
significantly between the two analyses with the monthly median amplitudes 
always smaller (e.g., 5.3 versus 8.8 gammas for the amplitude of the 
annual wave in H at College). Moreover, the corresponding statistical 
error estimates were always smaller in the case of monthly median data, 
indicating less interannual variability of the periodic waves when monthly 
medians are used. Thus, the following analyses are based on time series 
of monthly median values of H and Z. (Note that monthly mean MRN data 
were used here, as previously, because the range of those- fluctuations 
is relatively much smaller.) 

Monthly data of both the MRN and geomagnetic parameters were analyzed 
with the joint periodic regression technique of- Belmont, et al., 1974a. 

The technique can be used to analyze a time series of irregularly spaced 
data points, weighting the months by number of observations, even if zero, 
and to simultaneously determine an estimate of the statistical error of 
the amplitude and phase of each frequency included. Further, frequencies 
analyzed need not be integral divisions of the period of record. 

Frequencies included in the present analyses are the long-term mean, 11 
year cycle (geomagnetic data only), quasi -biennial oscillation (29 months), 
and the first four harmonics of the annual wave. 
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The statistical errors given in IJab le.. 2, .which .provide, cqnfidence . 
estimates for the results., are not the same as RIB -deviations . However, 
they resemble EMS deviations because .the coherence of the data from cycle 
to cycle is the most important consideration, ip,! computing .them;, .in- fact , 
regression of the errors in Table 2 (SE) for the annual and semiannual j. 
waves with RMS deviations determined from harmonic analysis of yearly 
data showed that SE = 0.38 RMS. ’The regression coefficient is'. small for 
Several reasons:' because the ' frequencies included' in- the SE -analysis -are 
not Orthogonal over the data, they *can interfere with each 'other'''t6'-give 
a better fit ' (-smaller residuals) to the complete time series than can^ 
the orthogonal components of harmonic analysis. Also, as SE 
weights each point of the time series by the number of .-i.o.bser;yations , it 
allows occasional erratic points based on few observations to be largely 
disregarded. 

2; Annual Wave 

Estimates of the -amplitude and phase, with errors, of .^the.,.Q^0, ^ 
annual, and semiannual waves in H and ,Z are given, in. Table ,2 . ,pResi|JLt,^^ for 
the annual .wave in H and Z are plotted in Figure 1 ,as functions ,of gftpuiag- 
netic. -latitude. Lines estimating the latitudinal variation in.^the- figure 
hav.e been fitted .by eye. These results are similar to..the corresponding 
values given by Currie (1966), but .have much less scatter', p.articularly, . 
at mid-latitudes. As noted by Currie, the scatter .of his, results ^may 
arise, from differing periods of .record at the various stations,, he .used' 
for example, the interannual variations of the annual..waves i-n;H. af Tqcson 
and Sitka (Figure 2) are so large that averaging a given number of .arbi- 
trary years -will clearly lead to widely varying mean values. In antici- 
pation of the discussion in Section B-, 'a' large part of the -interanriual 
variations in Figure .2 is due to the Well-known solar- cycle influence- i 
on E' region ionization. 'Note that both stations are at'^the- right extreme 
in ’1963' and 'the -left extreme in 1969 (1963 was near sun’ spot miiiiminn ’'and 
1969 near sun spot maximum). Returning to Figure 1, -the 'sharp x-riGreasej' 
of amplitude of the annual waves in H and Z at high latitudes has been 
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noted by Currie (1966) whose data extended to 80*^ (geomagnetic latitude). 
The decrease of amplitude in H and continued rise in Z as the pole is 
approached does not seem to have been reported previously, although 
Lange 1 and Brovm (1974) have noted that the largest seasonal variations 
of AZ are near the pole. 

The phase of the annual wave in H is fairly uniform at all latitudes, 
with the annual maximum occurring in June (Figure la). The phase of Z, on 
the other hand, undergoes an abrupt shift of 180° near 65°N. Equator- • 
ward of 65°N the average phases of Z and H are quite similar, 

3, Semiannual yJave 

Amplitudes and phases of the semiannual variations in H and Z are 
plotted in Figure 3. The present amplitude results are generally smaller 
than those given by Currie (1966). In this case, the difference may 
again be due to differing periods of record, or it may be due to our use 
of monthly median data which reduces the occasionally severe impact of 
magnetic storms which occur on a predominantly semiannual rhythm. It 
should be borne in mind that the present results are for the mean semi- 
annual wave over about one sunspot cycle. [Chapman and Bartels (1940) 
have shown that the semiannual amplitude varies with the sunspot cycle, 3 
From 70° to 80° magnetic latitude-, the decline of amplitude in H and the 
increase in Z are in accord with Currie's (1966) results which extend to 
Godhavn (80°N). As the pole is approached from 80°N the amplitudes of 
'both H and Z increase, although the large statistical errors associated 
with the H values make that analysis less reliable. 

The phase of the semiannual variation in H is fairly steady up to 
about 70°N, vith a 180° shift near 75°N indicated by all three stations 
north of 80°N. The phase of Z is less steady, but indicates a systematic 
shift with latitude such that the phase of the pole and the equator are 
about 180° different. 
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4. Quasi -Biennia I Oscillation (QBO) 

Before discussing the periodic analysis results for the QBO in H 
and Z, it should be pointed out that there have been several conflicting 
reports regarding the existence of a quasi -biennial line’ in the geomag- 
netic' spectrum. Hope (1963) reported that the QBO in K^’ had been 
isolated,’ but Currie (1966) could not find it in the spectra of H or Z 
and suggested' that these results were based on faulty numerical filtering 
procedures. Fraser-Smith (1972) presented the spectrum of the index 
and concluded that no QBO exists, but Currie (1973) has analyzed H and 
Z data from 49 observatories and now concludes that there is, a line near 
2 . 15 years . 

Nearly all periodic waves in geophysical data show variations from 
cycle to cycle, but usually the amplitude and phase converge on mean . 
values if enough cycles are averaged. Statistical tests pan be' used to 
determine if enough cycles of a periodic wave have been used to estimate 
the mean wave with confidence (Chapman and Bartels, 1940). As the quasi- 
biennial oscillation is not truly periodic, but has variable amplitude, 
phase, and period from cycle to cycle (e.g., see Figure 4), there is no 
assurance that a mean wave can be rigorously defined in a usual statis- 
tical sense. Thus, the mean QBO can only be defined for the years of 
record analyzed by each writer, with the understanding that the QBO for 
different years of record will probably not have the same amplitude or 
phase. For this reason, the latitudinal variation of the QBO values in 
Table 2 is erratic and inconclusive unless stations with the most complete 
and most nearly identical years of record are considered. Therefore, only 
those stations with over 110 months of data have been used in order to 
obtain the most reliable estimates of the latitudinal variation of the 
QBO. (The average period of the QBO during these years was 29 months.) 

I 

In Figure 5, the average QBO amplitude is near 2y for both elements, 
although Resolute (83°N) indicates an increase of H’s amplitude as the 
pole is approached. The individual phase, dates, relative to 1 January 
1960, have fairly uniform latitudinal variation except for Z at College 
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and H at Resolute. 

b; comparison of geomagnetic and meteorological periodic waves 

The present ohjective is to determine possible relations of geo- 
magnetic to meteorological variables by comparing periodic properties of 
geomagnetic and MRN data. Identifying those periodic frequencies which 
show a close relationship can allow effort to be focused on them, with 
the remainder of the variance discarded as being, unrelated and therefore 
of no immediate interest. 

1. Procedure 

Waves of the same period whose relative phase lags show broad 
patterns of spatial continuity are sometimes found to be related, and 
charts of the relative phase lags between MRN and geomagnetic periodic 
variations will be presented below. However, as all periodic components 
of the two data sets may have large year -to-year variability in ampli- 
tude and phase (for example, the annual wave in H in Figure 2), it is 
desirable to first examine the year-to-year relationship of each frequency 
to determine if " average— phase lag values are representative. The 
coherence square (C0H2) statistic of cross spectral analysis provides an 
objective measure of how uniformly the amplitudes and phases within each 
frequency band vary with time at a given location, and has been used 
here to decide which frequencies of MRN and geomagnetic data are sjmchro- 
nous,' Cross-spectral analyses of horizontal and vertical components of 
the geomagnetic field versus the temperature and wind at nearby rocket 
stations were made using monthly data at five stations with the most 
complete periods of record (Table 3), yet well distributed in latitude, 
and with a maximum lag of twelve months. Prior statistical significance 
of each C0H2 value was tested by the method of Julian (1975). No values 
in the frequency band centered at 24 months (near the QBO frequency) , nor 
more values than expected by chance for frequencies higher than 2tt/6 months, 
passed the 5% confidence level. C0H2 values for frequency bands centered 
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at biennial, annual, semiannual, and terannual periods are given in 
Table 3; those which exceed the 5%, 1%, and 0.17o confidence levels are 
marked . 

* c 

There is apparently closest coupling for the annual and semiannual 
variations of zonal wind with H at mid- latitudes and with Z at lower 
latitudes. The annual temperature variation, especially at 30 km, also 
is significantly coupled with geomagnetic variations . Semiannual varia- 
tions of temperature and all terannual variations exceed the 5 % confi- 
dence limit no more often than expected by chance. These results indicate 
that only the annual, and some semiannual, variations in MEN and geomag- 
netic data exhibit significantly synchronous year to year changes; . 
therefore, only those waves will be considered further. 

There are two techniques for determining relative phase lags: first, 

relative phase lag can be found during cross-spectral analysis. Second, 
the phases determined by periodic analyses can be subtracted. The latter 
method has the advantage that a measure of confidence can be derived by 
combining the statistical phase errors determined during periodic 
analysis. This was done by the root-sum-square technique. It was found 
that all phase lags associated with a C0H2 in Table 3 which exceeded the 
1% confidence limit were within the limits of statistical error of the 
phase lags determined by subtracting periodic analysis results. Therefore, 
values presented below are based on periodic analysis results. 

2 . Relative Phase Lags 

Relative phase lags for the annual and semiannual waves between MKN 
and geomagnetic data are presented in Figure 6 .as functions of height and 
latitude. For each station pair listed in Table la, the relative phase 
lag of each frequency for each parameter was determined by subtracting 
the phase of the geomagnetic wave from the meteorological wave, at 4 km 
height intervals from 28-64 km. The resulting phase lag values were 
plotted at the geographic latitude of the MRN station. Contours were 
drawn for phase lags = +30, +150 degrees to indicate areas of nearly in 
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or nearly out of phase. The relative uncertainty of each value, esti- 
mated by the root-sum-square of the individual phase errors, and the 
spatial patterns of phase give'n in Figures 1 and 3, and in Belmont, et al. 
(1974a), and Nastrom and Belmont (1975), were taken into account while 
drawing the contours . 

In Figure 6, U-H are out of phase throughout the mid-latitudes for 
both the annual and semiannual waves . U-Z are out of phase from about 
10°N - 40°N for the annual wave, with small in phase areas at high lati- 
tudes. The phase lags of TJ-Z for the semiannual wave are near 180° in 
the upper tropical stratosphere and nearly in phase at high latitudes. 

The annual waves in T-H are out of phase above 55 km near 20°N, and 
in phase north of a line from 28 km, 10°N to 64 km, 60°N. The annual 
waves in T-Z appear out of phase in the upper low-latitude stratosphere 
and at highest latitudes, and are in phase near 30-50°N. Clearly, the 
phase lags presented in Figure 6 have broad spatial continuity. Together 
with the large G0H2 values these results suggest that physical coupling 
between the MRN and geomagnetic periodic variations may exist. Possible 
mechanisms which could produce coupling will be discussed next. 

3. Discussion 


a. Annual Wave 

It must be noted here that the annual variation in geomagnetic 

data is not yet fully understood, although several writers have discussed 

it. Vestine (1954) suggested it could be a seasonal effect induced by 

air motions in the ionosphere. Currie (1966) concurred with Vestine and 

offered qualitative arguments from the scanty data then available, and 

later (Currie, 1974)- strengthened the theory by arguing the annual wave 

could not arise from modulation' of the S current system but must be a 

9 

DC effect. 

(1) Suggested Mechanism ; Due to the differing ion and electron 
Hall conductivities, zonal wind in the lower ionosphere produces a ring 
current along the wind which induces a magnetic field in the meridional 
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plane. This induced magnetic field affects the geomagnetic field in 
proportion to the wind speed, and the effect decreases with distance. 

The maximum effect on the N-S component of the geomagnetic field will be 
directly below and above the wind jet where the. induced field is, 
coincident with the geomagnetic H field. Similarly, there is a maximum 
effect on the Z component to the north and south of the zonal wind jet 
with a minimum directly below and above it. 

At high latitudes (Fig. 1), the maximum amplitude of -the annual 
wave in both H and Z occurs. It apparently has not yet been explained, 
in the literature. The cause .of the high latitude maximum could be the 
annual variation in ionization density, which is a function of solar 
elevation angle, and winds in the lower ionosphere or of magnetos pheric 
origin; but there is insufficient data to verify either hypothesis at 
this time. 


At mid- and low- latitudes , however, sufficient data are now 
available to crudely estimate the magnitude of the annual effect of iono- 
spheric winds on geomagnetism and thereby perhaps bring future research 
efforts on this issue into focus. Here the annual variation in wind is 
the major factor as there is only a small seasonal change in electron 
density. The annual -variation in zonal winds in the lower ionosphere 
has maximum amplitudes of about 30 m/s from 20-50° latitude near 110 km 
with phase dates near mid-May (Groves, 1972). Ionized gas is dragged 
eastward during the half year centered about May, and westward .during 
the half year centered about November, producing an annual variation in 
the geomagnetic field intensity. To estimate the magnitude of this 
effect, the current sheet approximation is applied using a width of 

500 km (after Bates, 1975), depth of 15 km, uniform charge density of 
4 -3 

5 X 10 cm , at an altitude of 110 km. A wind variation of (30 m/s) cos 
(wt+0) yields a field variation of (3.28 gammas) cos (wt+0)-. Of course 
this estimate could easily be changed by a factor of -two or more, but '• 
the amplitude is certainly of the proper order of magnitude for the mid- 
latitude annual wave in geomagnetism (Fig. 1), Also, the charge density 
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varies, particularly with solar zenith angle; this could account • 
windS' (mid-May) and the geomagnetic variation (mid -June) . Since the 
annual amplitude in the zonal wind above 100 km has a maximum near 50*^ 
it should create a maximum in the annual amplitude of H near 50*^ and a 
minimum in Z near 50%, as found in Eigure 1. 

Conventional heat sources (e.g., radiative heating) 'are adequate 
to account for the annual wind waves in the stratosphere (Leovy, 1964) 
and lower thermosphere (Volland and Mayr, 1972). The MRN individual 
data' has high coherence with the geomagnetic data at the' annual frequency 
because the variations in annual wave between the thermospheric and 
stratospheric wind are apparently also coherent. (The fact that the 
annual wave in the stratosphere is out of phase with that in the thermo- 
sphere has no bearing on their coherence.) The point here is that a 
seemingly intriguing relation between two parameters may arise from a 
mutual association with a third parameter through normally accepted 
processes; in this case the third parameter is the annual wave in thermo- 
spheric circulation. Hence, the present results do not suggest any 
geomagnetic influence on the atmospheric circulation. 

These results should be useful to those trying to understand 
apparent correlations between atmospheric and geomagnetic processes, and 
to those concerned with the description of the earth’s magnetic field 
and its variations . The annual dynamo concept presented above could be 
incorporated into models of -the geomagnetic field and thereby help 
overcome the problems of interpretation discussed by Alldredge and 
Stearns '(1974) . 

As the above calculation, based on constant ion density, does 
not pertain to the large, high latitude annual waves in geomagnetism, it 
is not inconsistent- that relatively low C0H2 values for the annual 
frequency are found in Table 3 at .Greely and Churchill. Values of G0H2 
■ in Table 3 at mid- and low- latitude stations are less than 1.0 for reasons 
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besides instrument error and incomplete sampling. Solar cycle 
influence on charge density in the ionosphere, causing a solar cycle in 
the annual wave in geomagnetism but not in that in stratospheric wind, 
may be the most important additional reason. However, upward propagating 
planetary and gravity waves, which may affect the stratosphere and iono- 
sphere much differently, could also be important. Although recent 
theories suggest that planetary waves will be absorbed, reflected, 
refracted and radiatively damped in the stratosphere and mesosphere, 
there is a large body of evidence which suggests they do exist in the 
lower thermosphere (e.g., Lysenko, et al., 1974; Deland and Friedman, 

1972; Graznik, et al., 1975). The possible role of gravity waves in the 
upper atmosphere is also poorly understood (Muller and Kingsley, 1974). 

It seems unlikely that these uncertainties will be cleared up until 
detailed wind measurements from the surface to the lower thermosphere 
are studied. A preliminary effort has been made by Manson, et al. 

(1975 ), but conclusive results are not yet available. 

(2) Applicability to Correlation Studies : During an early 

phase of the present study the linear correlation coefficients between 
the monthly means of MKN and geomagnetic data were computed. Those 
results, given in Table 4a, have a high level of statistical significance. 
It is now realized that the correlation coefficients are large because 
the annual waves in MRN and geomagnetic data are coupled and, except at 
low latitudes, the annual wave is generally larger than any other periodic 
component in the MRN data. Thus, one would expect the linear correlation 
between zonal winds and geomagnetic data to decrease significantly if 
the annual waves were removed from both data sets. To test this hypo- 
thesis, the linear correlation coefficients were recomputed between the 
monthly residuals after the annual waves had been subtracted. The results 
of this test, given in Table 4b, show that in nearly all cases the 
correlation ceases to be significant when the annual wave is removed. 

The correlation remains significant at Hawaii because the semiannual 
wave in zonal wind is nearly as large as the annual. 
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Application of this point to other reported correlations may 
help explain them. For example, King (1975) has reported that the 
longitudinal variations at 60°N of the average 500 mb height for 
January and the geomagnetic intensity shifted 25° in longitude have a 
correlation coefficient of -0.963. Longitudinal variations in the- 
circulation of the mid-stratosphere reveal a standing wave up to at 
least 10 mb'; in the meridional component the predominant standing 
wavenumber is two (van Loon, et al., 1972; Figure 72). Lysenko, et.al,. 
(1972), and Glass, et al. (1975), have offered evidence that standing 
waves also exist in the circulation of the lower thermosphere. If the 
predominant wavenvimber of wind speed, ion density, or a combination of 
them in the lower thermosphere in January is two, then the resulting 
current will induce a wavenumber two pattern in the longitudinal varia- 
tions of the geomagnetic field intensity. The high correlation found 
by King may therefore reflect a very mundane relationship, as long 
proposed by Wulf (1945), rather than any solar-terrestrial effect. A 
similar principle could apply regarding the relationship between spatial 
variations of tropospheric temperature, htomidity, and surface pressure 
and the geomagnetic field intensity reported by King (1974). 

b. Semiannual Wave 

The results in Table 3 and Figure 6 suggest that the semi- 
annual waves in MKN zonal wind and geomagnetism are also closely coupled. 
For the zonal wind in Table 3, significant C0H2 values are found for 
the semiannual variation at nearly the same station- levels as for the 
annual variation. A dynamo mechanism might be suggested, as Groves 
(1972) shows that there are large semiannual wind variations near 115 km. 
However, Volland and Mayr (1972) found that most of the latitudinally 
varying part of the semiannual wind wave in the lower thermosphere is 
due to corpuscular heating. They suggest (Mayr and Volland, 1971) that 
this heating is related to the semiannual occurrence of magnetic storms , 
which Chapman and Bartels (1940) have argued is due to earth-sun 
geometry and thus is independent of meteorological influence. There- 
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fore, the close coupling of the semiannual waves seen in Table 3 and 
Figure 6 may be explained independently of the dynamo mechanism of 
the annual wave. 

More insight regarding the cause for this coupling of the 
semiannual waves in geomagnetic and MRN data might be possible if the 
cause of the semiannual wave in zonal wind were known. Possible causes 
for the tropical semiannual wind wave have been discussed by Dickinson 
(1975), but the extra tropical semiannual wave has not yet been explained. 
As processes which show more symmetry in one coordinate system than 
another may be driven by mechanisms peculiar' to that coordinate system, 
tests of the relative symmetry of the semiannual wind wave in geomagnetic 
and geographic coordinates were made. These tests, described below, 
were generally inconclusive. Finally, three possible causes of the 
extratropical semiannual wind wave are discussed. - None can yet be 
accepted, and it is suggested that more research is needed before a 
conclusion can be reached. 

w 

(1) Further Tests for Coupling with the Geomagnetic Field : In 

order to test the relative symmetry of the semiannual wind wave in 
geomagnetic compared with geographic coordinates the relative phase 
lags of Figure 6b and 6d have been plotted in geomagnetic coordinates 
in Figure 7. In either case, the change of coordinates makes little 
difference, although for U-2 the contours become smoother in geomagnetic 
coordinates . 

Belmont, et al. (1974b), compared the symmetry of the ampli- 
tude of the semiannual wind wave at 50 km on maps in the geographic and 
geomagnetic coordinate systems and found the sjmimetry slightly greater 
in geomagnetic coordinates. Even greater symmetry may be found by plotting 
the amplitude at each station at that level where the closest relationship 
is found. The height of the level at each station was selected as that 
height where the magnitude of the product of the two semiannual waves' 
amplitudes and the cosine of their phase lag (a 2 *b 2 -00340) is maximum. 
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Note that this parameter, which is an approximation of the co-spectrum 
in the case of large C0H2, will be relatively- small if either amplitude 
is small or if the phases are near quadrature. The height of this 
surface is shown by the dotted line in Figure 6b. Contours of the ampli- 
tude of the semiannual wind wave at the heights thus selected are shown 
in Figure 8, and appear to show little, if any, enhanced sjraametry in 
either coordinate system compared with the results- of Belmont, et al. 
(1974b) . Clearly, these tests for increased symmetry do not suggest 
preference for either coordinate system. 

(2) Possible Mechanisms : Three hypotheses can be advanced to 

account for the extratropical stratospheric semiannual wave in zonal 
wind. Before discussing them, however , it should be pointed out that 
Gregory, et al. (1975b), have noted that the phases of the annual wind 
waves in the stratosphere and upper mesosphere are reversed. Cole and 
Kantor (1974) have noted a similar relationship with regard to the 
annual waves in temperature at stratospheric and mesospheric levels. - 
Both papers suggest that the semiannual waves in the lower mesosphere 
at extratropical latitudes result from the overlapping of the annual 
waves . This descriptive account of the lower mesospheric semiannual 
wave is useful, but does not by itself explain the semiannual wave. For 
example, early descriptions of the tropical semiannual wind wave in the 
upper stratosphere viewed it as the result of alternating intrusions 
of winter hemisphere westerlies into the summer hemisphere (Webb, 1966)-. 
While that does occur, it does not explain the tropical semiannual wave, 
and efforts to do so have invoked a wide variety of mechanisms, e.g., 
ozone heating, the diurnal tide, planetary waves, Kelvin waves, and 
semidiurnal tides. Hopefully, the discussion below will help stimulate 
other research efforts to explain the extratropical semiannual wind 
wave . 

The first forcing mechanism 'for the extratropical' semiannual wind 
wave to be considered here is upward propagating planetary waves. If 
these waves interact with the background flow oh a semiannual basis 
they could induce a semiannual component in the background wind speed. 
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In an effort to determine if the amount of absorption of planetary 
waves varies with season, the variance spectrum of filtered time series 
of Zonal and meridional winds have been determined at eight MRN stations 
on a seasonal basis. Details of procedure and complete results, are 
given in the Appendix. As noted in the Appendix, maximum power usually 
occurs between 2 tt/ 10 and 2 tt/ 20 days, so for brevity only the results at 
2n/ll days will be presented below. However, the graphs for the total 
variance and for the results at 2n/6.3 days (not shown) are similar. 

Planetary waves propagating vertically in a hydrostatic 
atmosphere with no dissipation increase their spectral density (i.e., 
power) exponentially with height. Further^ if attenuation of the waves 
occurs, the slope of the power will be proportional to the amount of 
attenuatipn. Spectral density for the frequency band centered at 2 tt/ 11 
days is presented as a function of height in Figure 9 for six MRN stations. 
Although the power at a given level changes with season, sometimes dramat- 
ically, the slope of the curves does not change much with season except 
at Kennedy. There is a large difference at Kennedy between the slope 
during spring and autumn from that during the solstitial seasons. During 
winter unexplained absorption occurs above 45 km at both Kennedy and 
Pt. Mugu. These results support the hypothesis of Belmont, et al. (1974b) 
that planetary wave absorption, is responsible for the secondary amplitude 
maximum of the semiannual wave found near 30°N. They do not, however, 
suggest that the semiannual wave at other latitudes arises directly from 
seasonal absorption of planetary waves. 

The second proposed mechanism is influence on the ozone field 
by particle precipitation. This mechanism was proposed by Belmont, et al. 
(1974b), but cannot yet be directly tested due to a dearth of high level 
ozone data. However-, it should be noted that Heath (1974) has found 
evidence for a non-photochemical source of high latitude ozone creation 
which he attributes to incident charged particles; and recent modeling 
efforts by Crutzen, et al. (1975), have shown that incident charged 
particles can dramatically influence the ozone field. Also, Golyshev, 
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et al. (1974)., found that the amplitude of the semiannual wind wave near 
the stratopause exhibits a solar cycle modulation. To illustrate this', 
yearly values of several solar and geophysical parameters are presented 
in Table 5. A' station-year is not included in the table unless data 
for all twelve months are available, and temperature data were thus too 
irregular to include in the table. Note that the values of the sunspot 
number and of the semiannual amplitudes have relative maxima in 1969 in 
all cases. Further, note that the annual wave in zonal wind is a rela- 
tive minimum in 1969 at all stations except at Barking Sands. While 
this table suffers from the short period of record available, it does 
support and extend the results of Golyshev, et al. (1974).^ 

Solar cycle modulation of the periodic variations in strato- 
spheric zonal wind, as seen in Table 5, is consistent with the hypothesis 
that particle precipitation during magnetic storms influences the ozone 
and hence thermal and wind fields . If stratospheric semiannual wind 
variations are related to the occurrence of magnetic storms through the 
ozone field, then their amplitude should be largest during active sun 
years (as found in Table 5) as the semiannual component in magnetic 
stoinn frequency is largest during active sun years. Solar cycle modula- 
tion of the annual wind wave is not easily conjectured, but the well- 
known solar cycles in total yearly magnetic storms and yearly solar flare 
occurrence may prove responsible, especially in view of the results of 
Crutzen, et al. (1975), regarding particle precipitation and ozone 
concentration . 

The third possible mechanism is IR radiation generated in the 
lower thermosphere during magnetic storms and absorbed by CO 2 at 30-40 km. 
During magnetic storms the amount of IR radiated by the lower thermosphere 
is increased by several orders of magnitude, and Gordiyets, et al. (1972), 
have suggested that it causes heating of CO 2 at 30-40 km and H 2 O at 7-12 
km. This mechanism has appeal because 30-40 km is the region where 
maximum amplitudes of the semiannual wave in observed temperatures occur 
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(Nastrom and Belmont, 1975), and the semiannual component in the occur- 
rence of magnetic storms would produce the proper phase and periodicity. 
Large amounts of radiative energy are possible for brief periods during 
severe geomagnetic disturbances, but following Volland and Mayr (1972), 
the long period form of this heat input (averaged over space and time) 
should take the same form as the variation in magnetic energy, which is 
given by 


U? ~ [1-0.2 cos(W ^t)], 

JL 


where U is a yearly mean magnetic energy, dependent on solar activity, 

W is the semiannual frequency, and t is time. This implies that the 

amount of energy deposited by IR radiation is E ~E [1-0,2 cos(W t)3, 

IK J.K S3. 

where E is a yearly mean value. Alternatively, because layer mean 
i.K 

temperature and energy are directly related, T ~T [1-0.2 cos(W t).]. 

‘ iJx IK S 3 

This latter relation says that the IR should contribute five times as 
much to the mean temperature as it does to the semiannual component of 
temperature. However, in order to produce a zonal wind oscillation of 
20 m/ s the latitudinal variation of the corresponding temperature 
oscillation must be near 5°K (Groves, 1972)-, which implies a contribu- 
tion to the average temperature of 25°K, As dynamic models of the 
stratosphere have encountered no evidence of such a large unconventional 
heat source (Leovy, 1964), it seems highly unlikely that IR radiation 
from the lower thermosphere is an important forcing mechanism for the 
stratospheric semiannual zonal wind wave. 

c. QBO 

Coupling between the QBO’s in MRN and geomagnetic data may exist 
despite the lack of statistical significance of the C0H2 values in Table 3. 
Even in the tropical stratosphere, where the QBO is the dominant oscillation, 
the QBO is not regular in amplitude or period from cycle to cycle nor 
between levels during the same cycle (Wallace, 1973). Thus, the small 
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C0H2 values may be misleading in this case. Moreover, the QBO in thermo- 
spheric zonal winds found by Sprenger, et al, (1975), suggests that the 
geomagnetic QBO may result from a dynamo mechanism, parallel to the 
annual wave. Although the present theory explaining the well-known 
tropical stratospheric QBO appears successful (Dickinson, 1975) it is 
dependent on waves and processes unique to the tropics, and thus- cannot 
be invoked to explain an extra tropical thermospheric QBO. Similarily, 
any explanation for the thermospheric QBO cannot be based on processes 
unique to the thermosphere because the large negative correlation between 
the multi-year variations in Z at Honolulu with the 56 km zonal wind at 
Barking Sands (Pig. 4) suggests the oscillation is not unique to thermo- 
spheric (dynamo) altitudes. Until the altitude and latitude progression 
of the QBO throughout the upper atmosphere is better known, no conclusion 
regarding the present results seems warranted. 

C. SPECTRAL CHANGES DURING SUDDEN WARMINGS 

In the stratosphere, rapid changes in the number, amplitude and 
phase of planetary waves are the major events during winter. These 
changes are sometimes associated with "sudden warmings" and it seems of 
interest to study the changes in MEN and geomagnetic parameters during 
these' disturbances . For this purpose , all years have been catagorized 
as either major sudden warming years (SW) or as other years (MSW) . A 
sudden warming is defined to occur when there is a "reversal of the 
polar circulation at 10 mb. (30 km) or below". During 1961-72, SW' were 
in 1962-63, 1965-66, 1967-68, 1969-70, and 1970-71, according to a list 
by Finger. 

Several different methods could be used to study the changes of 
parameters during SW. For example, as planetary wave activity and 
other events associated with a SW are global in nature (Quiroz, et al., 
1975) , spatial wavenumber analysis of global data may be used to detect 
changes in the planetary wave patterns. However, the MRN and geomagnetic 
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data are not sufficiently distributed geographically to permit detailed 
spatial analysis. Superposed epoch studies are often useful for single 
station analysis, but in the case of sudden" warmings it is difficult to 
meaningfully define a key -date. Indeed, the criteria used for defining 
the occurrence of a SW are admittedly arbitrary. Thus, the approach 
used here is to perform- power spectrum analysis of single station data 
and to compare the spectra of SW years with those of MSW years. 

The available wind observations (Hook, 1972; Gregory and Manson, 

1976) indicate that the circulation of the lower thermosphere is disturbed 
during a SW. Winds in the ionosphere can act as electric currents and 
can thereby produce variations in ,the geomagnetic field. Of course, 
processes unique to the magnetosphere can. also produce variations in the 
geomagnetic field; but if a geomagnetic spectral feature can be asso- 
iated with a meteorological process, it may be .reasonable to assume that 
it arises from that meteorological process. Thus, studying spectral - 
changes in the geomagnetic field between SW and MSW years may help 
better understand the thermosphere. Spectral analysis results for the 
zonal winds and for the horizontal field- intensity are presented first, 
with a brief discussion of .noteworthy features. A comparison of the two 
sets of results follows and a possible interpretation is suggested. 

1. Stratospheric Zonal Winds , 

Spectra for the zonal winds at 40 km at Fort Greely and White Sands 
are given in Figures 10-11 for autumn through spring. These stations 
were chosen because they have the most complete data at high- and mid- 
latitudes, respectively. In Figures 10-13, K is wavenumber, solid lines 
are 'for SW spectra, and dashed lines are for MSW spectra. In autumn and 
winter there is more energy at Fort Greely (Fig, 10) during SW years if 
the peaks near K = 8 are momentarily disregarded. The high frequency 
peaks will be discussed later. At White Sands (Fig. 11), however, 
largest energy occurs during MSW years at K '= 6 to 9 in autumn and K = 

2 to 6 in winter. In spring the high frequency energy is significantly 
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larger at both stations during years. Chi-squared confidence 
limits have been entered at noteworthly wavenumbers in the figures to 
indicate the probability that the differences arise from chance. 

These results for autumn and winter support Matsuno's (1971) suggestions 
that there is enhanced upward flux of wave energy at high latitudes 
(e.g., Greely) during SW years, but during MBW years the waves are 
refracted toward lower latitudes (e.g., WSMR) resulting in more energy 
there during IS^W years. 

2 . Horizontal Field Intensity 

Time spectra of the variations in H at College and Tucson are given 
in Figures 12 and 13. Spectra for several observatories were computed,. 
As College and Tucson demonstrate the salient features noted and are 
near the MRN stations used above, only they are presented here. The 
spectral differences between SW and MSW years noted below are probably 
due to meteorological influences, and not to solar induced effects. 
Hauska, et al. (1973), found that over all years 1932-1969 the magnitude 
of geomagnetic variations in the time range 4 to 40 days varies primarily 
with the approximately 11-year cycle. During the period 1961-1972, 

SW years defined above are well distributed over a solar cycle. 

In general, at both stations, there is either little difference 
between SW and >SW years or the spectral values are greater during MSW 
years. During autumn, the largest differences are at low wavenumbers 
(K=l~3 at CO and K=l-6 at TU) , while in spring differences are found at 
intermediate and high wavenumbers (K=4-5 at CO and K=3-ll at TU) . During 
MSW autumn at College, a significant (1% level) peak is found at K=9; 
less significant peaks at K=10 are found there during winter and spring. 

3. Discussion 

The following chart summarizes which years have the signifi- 
cantly greater spectral values and the wavenumbers at which they occur: 
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ATJT 

WIN 

SPR 

Greely 

SW(3-5) 

SW(5) 

MSW(3-11) 

CO 

MSW(l-3) 

MSW(l,5-6) 

None 

WSMR 

MSW(6-11) 

None 

MSW(6-9) 

TU 

MSW (1-6) 

None 

MSW(3-11) 


From the above chart and- Figures 10-13 two points should be noted. 
First, the only time the spectral values are significantly higher during 
SW years is at Greely during autumn and winter. The first point was 
noted to be consistent with the theory. of wave propagation and refraction. 
Second, significant spectral peaks near K=9 are found only at Greely and 
College during autumn and winter of MSW years- and at Greely during autumn 
of SW years. The second point may also be explained by planetary waves 
as will be suggested next. 

Planetary waves occur in the troposphere every year. As they propa- 
gate upward, they may be refracted toward lower latitudes or they may 
continue propagating upward, depending on the vertical and horizontal 
curvature of the flow profile (Simmons, 1974). As waves travel upward 
they decay; the rate of decay depends on the prevailing circulation. 

It is now hypothesized that waves of period near 4 to 5 days (K=8-ll) are 
upward propagating at high latitudes during all years . If during MSW years 
they do not suffer severe attenuation then they may continue all the 
way to the lower ionosphere resulting in spectral peaks near K=9 at 
Greely (40 km winds) and College (lower ionospheric winds). During SW 
years the prevailing circulation may cause large attenuation or total 
absorption; thus, during SW years the peak near K=9 at Greely in autumn 
is smaller than during MSW years , and a corresponding peak is not found 
in winter at Greely nor in autumn or winter at College. 

The above arguments, although sketchy and heuristic, are consistent 
with present knowledge. Lacking from present knowledge, however, is an 
adequate climatology of circulation differences during SW and MSW years, 
especially in the upper stratosphere. 
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D. 


CORRELATION OF MRN TEMPERATURES WITH K AND THE SOLAR SECTOR STRUCTURE 

P 


Numerous authors have suggested that the middle atmosphere may be 
heated following geomagnetic disturbances (e.g., Gordiyets, et al. , 1973). 
A desirable method of testing this hypothesis would be a superposed 
epoch study using a magnetic storm parameter as the keydate. However, 
as the MRN data are too scanty to permit that study, the alternative 
procedure of finding lagged correlations of with respect to the MEN 
temperature observations was used. The linear correlation coefficients 
between and the layer mean temperature, 40 - 50 km, at Fort Churchill 
are given in Figure 14. All temperature soundings, 1960-1972, which 
had data through the entire layer were used in this study. Values of K^, 
obtained from the World Data Center, Boulder, are' reported for three-hourly 
periods; thus, the correlation coefficient was determined at three-hourly 
intervals as the temperatures were lagged with respect to the values . 

In Figure 14, negative lag means that the value was measured before 
the temperature value. The relative maximum correlations are found at 
lag zero and at lag -15 hours; although both peaks are statistically 
significant at only the 5% level (if all data are assumed independent), 
these results are complementary to those found by others. 


Ramakrishna and Seshamani (1973) report a statistically significant 
correlation between the layer mean temperature (from 'grenade data) at 
Churchill, 60-89 km, and K^. The peak correlation occurs when tempera- 
ture is lagged 15 hours, and the largest correlation coefficients are 
found when the mean temperature is for the entire layer rather than 
just the upper portion of the layer. They report this correlation is 
significant at the 0.17o level. The largest regression coefficients, a 
measure of the relative magnitude of the effect, are found when only 
the upper portion of the altitude layer is used; and they suggest this 
may indicate a larger heating effect at highest altitudes. However, 
it also may be due to cancellation resulting from opposite effects in 
different portions of the entire layer. Results given by Zadvernyuk 
(1973) indicate that the critical' layers of the atmosphere may respond 
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differently to magnetic disturbances; e.g,, following a magnetic storm 
there may be heating at the mesopause but cooling at the stratopause. 

Several possible mechanisms could be suggested to account for the 
correlations discussed above: e.g., corpuscular heating, enhanced IR 

radiative exchange, etc. However, the correlation could also arise 
from a meteorological influence on K in a manner similar to that suggested 
by Hines (1973). Therefore, if correlation studies such as the above 
are to be taken as indicators of a geophysical process important to the 
lower atmosphere, they must be based on unambiguous parameters so that 
cause and effect can be clearly discerned. 

The previous studies are inconclusive with respect to solar-terres- 
trial effects. Wilcox (1975) has already related K to solar sectors. 

The real question is whether temperature can be related to solar sector 
structure. To examine this, a superposed epoch study of the 40 - 50 km 
layer mean temperatures with solar sector boundary crossings used as 
key dates is desirable. 

However, well defined solar sector boundaries sweep past the Earth 
at irregular intervals, about every week on the average, and the joint 
distribution of them with the intermittent MRN observations is not 
adequate for a superposed epoch study. Thus, it is possible to present 
in Figure 15 the sign of the temperature change at 40 km between 
closely spaced consecutive MEN observations at Churchill as a function of 
time, relative to a solar sector boundary crossing and K . The magnitudes 
of the temperature changes are not shown, but they are random. Solar 
sector boundary crossing dates were taken from the list in Shapley, et al. 
(1975). There are 21 temperature rises and 28 temperature falls on the 
chart, and they seem to be evenly distributed on both sides of the 
boundary. From this small sample it appears that the temperature trend 
shows no preference relative to the passage of a solar sector boundary. 

In summary, either the purported T-K^ relationship is due to some factor 
other than mutual coupling with solar sector structure, or a much larger 
sample would be required to establish reliably such a relationship. 
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IV. SUMMARY 


Periodic analysis results of the horizontal and vertical field 
intensity show that maximum amplitudes of semiannual and annual waves 
are at high latitudes. It is suggested the high latitude maxima of the 
annual waves arise from annual waves in ionization density and thermo- 
spheric zonal wind speed. At mid- and low-latitudes the annual' wave in 
zonal wind speed in the lower thermosphere (the dynamo region) , which 
is driven by solar heating, accounts for most of the geomagnetic annual 
wave . 


Annual variations in geomagnetic and MRN data are closely coupled. 

In view of the above cause of the geomagnetic annual wave, coupling 
between the circulations of the stratosphere and lower thermosphere can 
explain the geomagnetic-MRN coupling; thus, the present results for the 
annual wave are not evidence of any geomagnetic influence on the lower 
atmosphere. Other apparent correlations of the lower atmosphere and 
geomagnetic field may arise from a similar dynamo action in the thermo- 
sphere caused by coupling of the th ermosphere with the lower atmosphere. 

Semiannual variations in geomagnetic and stratospheric zonal wind 

data are also closely coupled. As the semiannual wave in thermospheric 

zonal wind is driven primarily by auroral heating, the cause of this 

coupling is not clear. It would be helpful if the cause of the strato- 

sphereic semiannual wind wave were known, so three possible causes were 

discussed. Planetary wave absorption seems to be a direct cause only 
o 

near 30 and heating by IR from the lower thermosphere during magnetic 
storms is energetically unlikely. Possible modulation of the ozone (and 
hence thermal and wind fields) by particle precipitation during geomag- 
netic storms has not yet been verified by observations . Amplitudes of the 
annual and semiannual waves in stratospheric zonal wind may be modulated 
with the solar cycle as they generally have extreme values concurrent 
with extreme values of the sunspot number. This result for the annual 
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■wave is believed reasonable as the yearly number of proton solar flares 
varies with the solar cycle and solar flares can dramatically affect the 
ozone field. If particle precipitation during geomagnetic storms also 
influences the ozone field, then this result for the semiannual wave could 
also be explained as the semiannual variation in geomagnetic activity 
varies with the solar cycle. 

Power spectrum analysis of zonal wind variations shows that at 
high latitudes there is significantly more wave energy in the upper 
stratosphere during years when major sudden warmings (SW) occur, but at 
raid- latitudes largest wave energy is found during years when major sudden 
warmings do not occur (MSW) . This could be explained by wave refraction 
which occurs in varying degree each year depending on the profile of the 
background flow; however, a climatblogy of background flows during SW 
and MSW years is apparently not available. If geomagnetic variations 
reflect wind activity in the lower themosphere then the noted differences 
between SW and MSW years at mid- and high-latitudes seem consistent with 
recent theories of planetary wave propagation. The present results 
suggest that the planetary wave absorption peculiar to SW years occurs 
in the upper stratosphere, far below the region where direct geomagnetic 
effects are significant, and thus any direct geomagnetic "trigger" for 
sudden warmings seems unlikely. 


Although the correlation between stratospheric temperature and K 

P 

appears statistically significant and is complimentary to the results of 
others, cause and effect cannot be discerned. If correlation studies are 
to be used as evidence of a solar-terres trial effect, they must be based 
on parameters of strictly solar origin such as the solar sector structure. 



28 


APPENDIX . 


Organized wave activity in the upper stratosphere has been studied 
with MRN data by several writers, most recently by Hirota (1975). The 
latter used only those MEN stations which had at least 30 observations 
during a given season, subjectively interpolated the data to daily values 
by analyzing height-time sections for each station, and computed the frequency 
content of the interpolated data by power spectral analysis. Hirota* s method 
is very effective for analyzing a single season’s data; however, the present 
objective, is to, prepare a climatology of the power spectra of MRN data, and 
a less restrictive, objective approach is desirable. Rocket data have 
historically been taken on an irregular often sporadic basis, and there 
are instances of many observations at a given station over a time span of 
a few weeks with relatively sparse, data before and after that period, A 
climatological analysis method should take advantage of those intermittent 
periods of dense data. The lag-weighted autovariance function method 
described below is suited for this purpose, and has been used to estimate 
seasonal power spectra of MRU wind components, 30-60 km. This method was 
also used to analyze the power spectra of geomagnetic variations reported in 
Section III-C of the text. Although a complete description of this method 
can be found in Dartt and Hovland (1974), a basic outline of it and the 
variations used here will be, given, 

A. BASIC DATA HANDLING AND TECHNIQUE 

At each 2 km level, 30 - 62 km, multiple rocket ascents over -a two-day 
period were averaged together and counted as one datiim in the time series. 

Due to the poor distribution of MRN observations, many data points represent 
only one observation and many are missing; but a surprisingly large number 
(for example, 20% at Churchill) of data points do represent multiple obser- 
vations. Interpolation was not used for miss.ing data. The time series thus 
obtained at each station and level were then high-pass filtered by convolution 
with a discrete, symmetric series of Gauss'ian .weights . To account for missing 
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observations, the weights under the filter were normalized at each data 
point such that their sum was always equal to 1.0. The ideal frequency 
response of this filter is shown in Figure A-1; -however, due to missing 
observations, the actual frequency response is slightly less sharp than 
shoxm in. the figure. 

If an observation is far removed in time from other elements of the 
time series the filtering process will be ineffective as the datum is then 
filtered, essentially, with only itself. To preclude this, it was required 
that there be at least five other data points under the filter (out of a 
possible 30) and that the sum of the weights before normalization be at 
least 0.25. These latter conditions resulted in discarding about 107® of 
the data. 

Autovariances up to Tag 11 were computed for each individual season, 
1961-1972. Three-month seasons were used at all stations with winter 
defined as December through February. The autovariances and the number of 
data pairs at each lag and each season were then stored for future use. 
Seasonal autovariances for all years of record were computed by combining 
individual seasonal values according to: 


_ 2 R(T)li(T ) 

R(t) S N(t) 


where R(t) is the autovariance and N(t) is the number of data pairs avail- 
able at lag T for a given season. With this procedure seasonal autovari- 
ances can be computed for the entire period of record, or for just selected 
years (e.g., years of major sudden warmings). Autovariances thus obtained 
were Banned J estimates of the power spectra were obtained by taking the 
cosine transform of the Banned autovariances . Finally, the computed values, 
V, were normalized: 

A 

S = V [2 MAX AT/2 it] 
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where, in this case, MA.X is 11 and AT is 2 days; and S is the normalized 
value. 

The percentage relative error of each spectral estimate was 
computed by determining the variance of Banned spectral estimates accord- 
ing to the formulation of Eddy (1968). The effective number of degrees 
of freedom required for that calculation were determined with the method 
of Mitchell (1963). It must be noted that these errors reflect how well 
each spectrum conforms to a particular statistical model and are 
only as realistic as that model. Further, they do not account for 
suspected error sources such as aliasing. Aliasing, or spectral folding, 
.results from sampling at a frequency lower than twice that of the natural 
variability; this problem is discussed in detail by Lumley and Panofsky 
(1964). If the true spectrum is a "red-noise" spectrum, as frequency 
decreases energy increases, then aliasing will tend to make the 
estimated spectrum flat, i.e., with equal energy at all frequencies. As 
discussed below,. this problem may be more serious for the meridional wind 
than for the zonal wind. 

B. TABLES OF SPECTRAL ESTIMATES 

Eight MRN stations have adequate data to provide meaningful estimates 
of the variance spectrum for the wind components. Temperature observa- 
tions are less plentiful than wind observations and did not provide use- 
ful results. Tables associated with this appendix give climatological 
spectral estimates for the wind components for each season. The values 
in the tables have been smoothed with height by a three point binomial 
filter. As the effect of missing observations on the frequency response 
of the high pass filter is difficult to estimate, no attempt to restore 
the spectra has been made; but it appears that the energy in the first 
frequency (centered at 2n/44 days) is reduced more than the 55% predicted 
by the theoretical frequency response of the filter. 

The relative reliability of the spectrum at each height is indicated 
by the number of lagged data pairs and the percentage relative error 
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estimates. The distribution of MEIN observations is such that N(l) . . . . 
N(ll) ~ 1/2 N(0), so only N(l) is given in the interest of brevity. Also, 
the percentage relative errors are nearly linear with wavenumber; thus, 
errors at intermediate frequencies can be estimated by linear interpola- 
tion of the values given for K=1 and K=ll. 

In spectral analysis there is always a trade-off between reaolution 
and reliability of the spectrum. By averaging adjacent bands a more 
reliable spectrum may be obtained, but a corresponding loss of resolution 
results. The most reliable parameter is thus the total variance of the 
filtered data, which is included in the tables. From experience, the 
best indicator of the reliability of the total variance is the number of 
lagged data pairs, and when N(l) is less than about 60 the variance 
should be disregarded. 

In the tables, "VAR" is the total variance of the filtered data, 

"N" is the number of data pairs at lag one, and "P.R.E," is the percent 
relative error for bands one and eleven. 

C. DISCUSSION 

During the Northern Hemisphere summer, the power of zonal wind 
spectra at a given level generally decreases with increasing K at mid- 
latitude stations (Kennedy through Wallops) and at Ascension. Zonal 
wind spectra at Greely, Churchill, and Barking Sands, and meridional 
wind spectra at all stations, are generally very flat at a given height 
but have large gradients with height. Although this strongly suggests 
that aliasing may be a serious problem at the latter stations, Dartt 
and Hovland (1974) report that summer spectra in the lower stratosphere 
(30 mb), detemined from relatively complete time series of twice daily 
radiosonde data, are also very flat, especially for the meridional wind. 

It is therefore likely that aliasing by periods longer than diurnal is 
not serious in the present results. Note that it is impossible to 
comment on possible aliasing by periods shorter than diurnal. 
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Figure 1, Annual wave in (a) horizontal, (b) vertical 
field Intensity, Dots are amplitudes and crosses are 
phases of stations in Table 2, plotted at geomagnetic 
latitude. Lines estimating the latitudinal variation 
of amplitude (solid) and phase (dotted) have been 
fitted by eye. 







Figure 2, Harmonic dial of the amplitude and phase of 
the annual wave for each year, 1960-1972, at Sitka and 
Tucson. The average values for all years are circled. 
The small number above each point is the year that the 
point represents (i.e., 3 is for 1963). Axes are 
labeled in gammas. 



Figure 3. As in Figure 1 except for the semi- 
annual wave. 




Figure 4, Twelve-month running mean values of Z (solid line) at Honolulu 
relative to the parabolic secular trend line, and twelve-month running 
mean values of zonal wind at 56 km at Barking Sands (dotted line) . Tick 
marks on abcissa are for 1 July; every third month is plotted. 
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Figure 6b. 























Figures 8a and 8b, The amplitude (m s~^) of the semiannual 
wave in zonal wind at the altitude indicated by the dashed 
line in Figure 6(b) or Table 1. Dots are MRN station loca- 
tions. (a) geographic coordinates; (b) geomagnetic 
coordinates . 
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Figure 10. Power spectrum of zonal 
wind variations at 40 km at Fort Greely’. 
Solid line is for SW years, dotted line 
is for MS¥ years, as defined in the text, 
Statistical significance level of the 
difference between the two spectra is 
indicated for continuum values . 

(a) autumn (b) winter (c) spring. 


Figure 11. As in Figure 
10 except at White Sands. 
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Figure 14. Linear correlation coefficient 
between layer mean .temperature, 40-50 km, at 
Fort Churchill and K as a function of the 
lag of temperature. 



Figure A-1. Theoretical frequency response of 
the numerical filter described in the Appendix. 
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Figure 15. Trend of the temperature at 40 km at Fort 
Churchill between closely spaced ascents made near a 
solar sector boundary crossing. End points of each 
line segment are plotted at the time relative to 
boundary crossing and at the corresponding value 

of K . 

P 





Table 1. List of meteorological rocket -stations . 


STATION 

"k 

LAT. 

k 

LONG. 

* 

LAT. 

LONG. 

NUMBER OF 

OES 

. at ; 


(GEOGRAPHIC) 

(GEOMAGNETIC) 

(WIND) 


(TEMP, 

a. Stations used 

in Figures 6 

-8. 





THULE 

77 

69 

88 

10 

335 


296 

FORT GREELY 

64 

146 

64 

261 ' 

1011 


563 

CHURCHILL 

59 

94 

68 

324 

991 


884 

PRIMROSE UKE 

55 

110 

62 

305 

316 


312 

WALLOPS 

38 

75 

48 

351 

1351 


674 

POINT MUGU 

34 

119 

41 

302 

1971 


1226 

WHITE SANDS 

32 

106 

42 

317 

2481 


988 

KENNEDY 

28 

81 

38 

347 

1916 


1142 

BARKING SANDS 

22 

160 

21 

265 

1372 


898 

ANTIGUA - . 

17 

. 62. 

28 

10 

466 


371 

SHERMAN 

9 

80 

20 

350 

63 i 


422 

kwAJALEIH 

? 

-168 

1 

238 

318, 


305 

ASCENSION 

-8 

14 

.-1 

55 

1196 

, 

9.37 

b. Stations used 

only 

in Figure 8. 





HEISS ISIAND 

81 

-58 

72 

156 

156 

(56 

KM) 

WEST GEIHNISH 

57 

, 7 

1,60 . ■ 

■-,84 

124 

(56' 

KM).. 

VOLGOGRAD” 

49 ■ 

-44 

‘ - L 

43 

125 

87’ 

(52 

KM) 

RYORI 

39 

-142 

29 

207 

32 

(48 

KM) 

ARENOSILLO 

37 

7 

41 

j76 

80 

'(54 

KM) 

SONMIANI 

25 

-67 

16 

137 

54 

t « 

(50 

KM) 

GRAND TURK 

21 

71 

32 

357 

1 

170 

(50 

KM) 

THUMBA 

9 

-77,/ 

0 „ 

i46 

145 

,(50 

KM) 

natal 

-6,-"' 

- '-35 

' 5 

|34 

131 

(46 

KM) 


Minus is south or east. 


NEAREST GEOMAGNETIC 
OBSERVATORY 

SYMBOL 

RESOLUTE 

RB 

COLLEGE 

CO 

CHURCHILL 

FC 

SITKA 

SI 

FREDRICKS BURG 

FR 

BOULDER 

BD 

TUCSON 

TU 

DALLAS 

DS 

HONOLULU 

HO 

SAN JUAN 

SJ 

FUQUENE 

FQ 

. GUAM 

GU 

HUANCAYO 

HU 



TABLE 2, Periodic analysis results of the geomagnetic field elements. Amplitudes are in tenths of gammas and 
phases are in degrees. Statistical errors are in parentheses. 



GEO- 

MAGNETIC 

LAT LON 



NUMBER 

OF 

. MONTHS 
H Z 

QBO 

AMP PHASE 

H 

ANNUAL 
AMP PHASE 

SEMIANNUAL 
AMP PHASE 

QBO 

AMP PHASE 

Z 

ANNUAL 

AMP PHASE 

SEMIANNUAL 
AMP PHASE 

THULE 

89 

1 

358 

34 

34 

46(25) 

-78(44) 

47(17) 

-177(23) 

53(17)-175(20) 

190(69)-102(36) 

126(42)-151(20) 

138(40) 

151(18) 

ALERT 

86 

168 

77 

77 

46(16) 

106(22) 

36(15) 

21(27) 

39(15)-154(25) 

15(24) 84(9) 

156(23) 4(8) 

89(22) 

-173(15) 

RESOLUTE 

83 

168 

115 

115 

. 40(5) 

159(7) 

29(5) 

-62(10) 

13(5) 153(23) 

29 (13) -118 (30) 

192(14) 1(4) 

87(13) 

-125(9) 

BAKER LAKE ' 

74 

315 

101. 

101 

89(13) 

-15(9) 

169(14) 

174(5) 

32(13) 55(27) 

109(20) 134(11) 

112(20) 2(11) 

79(20) 

-94(15) 

LIERVOGUR 

70 

71 

87 

81 

1(13) 

-122(95) 

12(14) 

27(23) 

50(4) 20(is) 

18(14) -7(14) 

21(14) 23(11) 

17(4) 

-80(14) 

CHURCHILL 

69 

323 

69 

69 

86(14) 

79(8) 

102(12) 

-151(6) 

47(11) -8(14) 

28(12) 81(26) 

97(11) 54(7) 

9(8) 

141(73) 

BARROW 

69 

241 

95 

94 

29(9) 

-107(17) 

94(8) 

-167(5) 

64(8) 21(7) 

16(6) 150(25) 

36(6) 142(10) 

46(6) 

-141(8) 

GREAT WHALE R1 

V 67 

347 

64 

64 

139(16) 

30(7) 

51(5) 

-74(19) 

39(15) -33(25) 

40(9) -43(13) 

41(8) -135(12) 

73(9) 

-66(7) 

COLLEGE 

65 

257 

156 

156 

22(5) 

-43(15) 

53(6) 

164(6) 

58(6) 38(5) 

15(3) 112(13) 

19(3) 119(10) 

14(3) 

-40(14) 

LERWICK 

63 

89 

20 

20 

24(10) 

81(17) 

21(16) 

-164(26) 

31(5) 10(10) 

13(16) 87(20) 

17(4) 7(12) 

15(3) 

-68(12) 

MEAKOOK 

62 

301 

68 

68 

7(4) 

30(39) 

44(4) 

-168(5) 

37(4) 25(6) 

14(4) -35(15) 

15(4) 6(14) 

9(3) 

-115(25) 

SITKA 

60 

275 

156 

156 

15(4) 

-49(15) 

48(4) 

168(4) 

41(4) 30(5) 

25(4) -137(5) 

7(2) 21(21) 

9(2) 

-49(16) 

FREDRICKS BURG 

50 

350 

156 

156 

22(5) 

-70(15) 

64(5) 

163(5) 

34(5) 37(9) 

26(5) -148(13) 

3(4) -115(86) 

13(5) 

-69(25) 

BOULDER 

49 

317 

72 

72 

19(5) 

-11(17) 

54(5) 

-174(6) 

33(5) 27(9) 

17(3) -140(12) 

12(3) 172(16) 

4(3) 

-80(61) 

STEPANORKA 

44 

111 

32 

32 

9(8) 

143(70) 

40(9) 

-175(14) 

32(9) -2(16) 

47(6) 104(9) 

8(5) 78(55) 

2(4) 

93(87) 

CASTLE ROCK 

43 

299 

33 

33 

- 19(6) 

-88(19) 

34(6) 

-170(11) 

36(6) -18(10) 

9(3) 101(19) 

6(3) -141(33) 

15(3) 

108(11) 

DALLAS 

43 

328 

99 

99 

15(5) 

-94(20) 

51(5) 

-168(5) 

23(5) 29(12) 

28(3) -96(6) 

17(3) 167(10) 

10(3) 

-62(18) 

TUCSON 

40 

312 

156 

153 

15(4) 

-39(18) 

53(4) 

152(5) 

32(4) 44(8) 

21(3) -158(8) 

13(3) -114(4) 

5(3) 

-43(44) 

SAN JOAN 

30 

3 

156. 

156 

9(6) 

-7(69) 

38(4) 

175(7) 

18(4) 43(14) 

26(4) -108(9) 

40(7) -138(10) 

28(7) 

-37(15) 

HONOLULU 

21 

266 


144 

15(4) 

49(15) 

35(4) 

143(6) 

19(4) 37(12) 

12(2) -116(9) 

25(2) -169(4) 

14(2) 

-7(8) 

FUQUENE 

17 

355 


38 

8(5) 

126(42) 

5(4) 

167(66) 

26(6) 85(12) 

81(18) 52(13) 

57(16) -83(37) 

53(26) 

-14(28) 

GUAM 

4 

213 

154. 

147 

15(4) 

-45(18) 

29(4) 

138(9) 

27(4) 84(10) 

10(3) -75(16) 

46(3) -174(4) 

11(3) 

-9(15) 

MUNTINLUPA 

3 

190 

60 

60 

84(7) 

50(5) 

18(7) 

-172(23) 

20(7) -31(21) 

19(4) -48(13) 

32(4) -162(7) 

22(4) 

-88(11) 

HUANCAYO 

-1 

354 

51 

44 

38(6). 

69(10) 

24(6) 

-138(16) 

4(4) 149(77) 

4(3) -13(58) 

33(3) 176(6) 

7(3) 

13(38) 

TOOLANGI 

-47 

221 

57 

0 

3(5) 

136(90) 

46(7) 

24(8) 

34(7) 50(11) 


- 

- 

- 

ARGENTINE IS. 

-54 

3 

96 

96 

10(3) 

52(22) 

36(3) 

32(5) 

33(3) 50(6) 

10(3) 13(18) 

26(3) 20(7) 

21(3) 

65(8) 

KERGUELEN 

-57 

128 

55 

0 

mSSi 

-19(16) 

26(6) 

2(14) 

34(6) 17(11) 

- 

- 

- 

- 

BYRD 

-71 

336 

91 

90 


62(5) 

115(35) 

95(19) 

15(25)-100(89) 

21(21) 61(5) 

73(19) 72(16) 

64(19) 

2(19) 











lable 3. Coherence-squared statistics between MRN and geomagnetic monthly data, and number of data pairs at 
lag one month, Station symbols as in Table 1, 


30 KM 48 KM 56 KM 


STATION 

CO 

FC 

FR 

TU 

HO 

CO 

FC 

FR 

TU 

110 

CO 

FC 

' FR 

TU 

HO 

N(T) 

95" 

65 

100 

134 

87 

86 

65 

89 , 

130 , 

85 

77 

65 

74 

114 

77 

N(U) 

114 

66 

144 

148 

115' 

110 

66 

144 

147 

117 

102 

61 

131 

146 

110 

a. H-T 

.12 

.50 

.23 

.50* 

.12 

'.Zl 

.54 

.15 

.19 

,02 

.16 

.56 

,16 

,20 

.08 

Z-T 

.05 

•51 

.06 

.18 

.46 

S3S 

.14 

.55 

.12 

.05 

.08 

.09 

•15 

,09 

.10 

.31 

H-U 

.16 

.32 

.48* 

.52* 

.27 

.20 

.59 

.51* 

.49* 

.27 

.15 

.53 

.48* 

.51* 

.28 

Z-U 

.03 

.37 

.00 

.13 

.53 

.02. 

•53 

.00 

.15 

.57* 

.02 

.55 

.00 

.16 

.56* 

b. H-I 

.39 

.28 

.05 

.11 

.13 

.19 

.18 

.12 

.05 

.18 

.11 

.10 

.16 

.11 

.14 

Z-I 

.18 

.05 

.19 

.17 

.34 

.32 

.18 

.17 

.03 

.08 

.05 

.05 

.36 

.13 

.02 

H-U 

.05 

.42 

•42 

•42 

.09 

.20 

.29 


.50>* 

.22 

.21 

.23 

■M 

.57* 

.27 

Z-U 

.15 

.25 

.17 

.07 

.48 

,16 

.03 

.07 

.04 

•53 

.20 

.10 


.03 

.66* 

c.. H-I 

.19 

.38 

.20 

.00 

.02 

.27 

.22 

.05 

.17 

.08 

.07 

.14 

.07 

.13 

.02 

Z-T 

.01 

.03 

.03 

.00 

.17 

.04 

.03 

.03 

.17 

.02 

.03 

.05 

.15 

.08 

.06 

H-U 

.01 

.04 

.03 

.12 

.31 

.01 

.09 

.03 

.07 

.05 

.02 

.03 

.02 

.07 

.09 

Z-U 

.20' 

.23 

.02 

.01 

.02 

.02 

.13 

.03 

.04 

.11 

.01 

.04 

.06 

.10 

.24 

d. H-T 

.02 

.28 

. .11 

.20 

.06 

.06 

.36 

.02 

.18 

.01 

,14 

.24 

.01 

.13 

,20 

. Z-T 

.22 

.36 

.03 

.07, 

.07 

.13 

.45 

.05 

.10 

.08 

.01 

.38 

.08 

.10 

.14 

H-U 

.05 

.12 

.08 

.13 

.16 

.08 

.32 

.09 

.11 

.03 

.06 

.37 

.09 

.10 

.04 

Z-U 

.05 

.38 

.02 

.03 

.07 

.02 

.35 

.00 

,.03 

.27 

.01 

.41 

.01 

.03 

.15 


CODE: a. Annual Wave 

b. Semiannual Wave 
,.c. 1 .Terannual Wave 
d. Quasi-biennial Wave 


STATISTICAL SIGNIFICANCE INDICATORS: * Value exceeds 0.1% confidence level. 

_ Value exceeds 1% confidence level. 

_ Value exceeds 5% confidence level. '• 


Station codes are given In Table 1. 

"H-T"' = Comparison of horizontal component of the geomagnetic field strength with temperature at 
• . specified ^ level from, nearest rocket, observation. 
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Table 4a. Linear correlation coefficients of monthly values of MRN and geo- 
magnetic data. Those which meet the 1% significance level are 
asterisked, 5% level are underlined. 


PARAMETERS 

(LEVEL) 

CO 

FC 

FR 

TU, ■ 

HO 

T-H 

30 

.083 

.-491* 

.223 

’ .280* 

-.138 


48 

.121 

,485* 

.014 

.247* 

-.070 


56 

.145 

.363* 

-.. 127 

.161 

.071 

U-H 

30 

-.188 

-.354* 

-.313* 

-.218* 

-.211 


48 

■ -.277* 

-.488*- 

-.350* 

^.320* : 

-.265* 


56 

-.338* 

-• -.466* 

-.343* 

-.351* 

-.222 

T-Z 

30 

.230 

-.280 

.117 

-.026 .. 

..293* 


48 

.253 

■ -.301 

.089 

-.181 

-.143 


56 

.118 

-.309 

-.164 

-.246 

-.303* 

U-Z 

30 

-.020 

.286 

.064 

-.110 ■ 

-.282* 


48 

-.026 

.284 

.047 

-.103 

-.490* 


56 

.010 

.303 

.051 

-.094 

-.504* 


Table 4b. 

Same as above 

except the annual waves 

were 

first subtracted from 


data. Note 

with above. 

the lack 

of significant 

correlation here compared 

T-H 

30 

-.114 

.174 

.111 

.099 

-.227 


48 • 

-.195 

.136 

.034 

.165 

-.044 


56 

-.024 

,006 

.090 

.195 

.173 

U-H 

30 

.013 

-.033 

.077 

-.038- - 

-.077 


48 

-.149 

-.226 

.083 

-.105 

-.051 


56 

-.190 

-.200 

.020 

-.109- 

.053 

T-Z 

30 

.169 

-.002 

.132 

-.024 

' .008 


48 

.154 

-.048 

.111 

-.167 

-.102 


56 

-.027 

-.151 

.015 

-.184 

-.172 

U-Z 

30 

.095 . 

.101 

.116 

-.026’ 

-.189 


48 

.083 

.086 

.020 

-.063 

-.203 


56 

.113 

.157 

.014 

-.082 

-.196 



TABLE 5. Yearly values 
wave relative 

o£_ solar and geophysical parameters 
minima are underlined. 

Some 

: relative maxim are 

underlined; 

for annual 

YEARS 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

SUN SPOT NO. 

• 54 . 

38 ' 

28 

10 

15 

47 

■94 

106 

106 ‘ 

105 

67 

10.7 cm FLUX 

104 

84 

80 

72 

76 

103 

143 

149 

151 ' 

156 

113 

SEMIANNUAL AMPLITUDES 












SITKA (H) 

4.7 •' 

4.3 

4.0 

3.2 

2.3 

4.1 

2.7 

3.7 

6'. 9 

5‘.S 

4.1 

FREDRICKSBURG (H) 

2.9 

3.9 

5.0 

4.1 

i.6 

4.5 

2.8 

i.-9 

6.9 

5; 5 

3.4 

TUCSON (H) 

3.3 

3.1 

4.4 

4.0 , 

2.3 

4.5 

4.8 

2.6 

6.1 

4.5 

2.6 

GREELY <U - 48KM) 





. 9.2 

9.2 

12.8 

19.5 

21.5 

11.9 


WALLOPS (U - A8KM) 

23.7 

“ 

24.0 • 


9.9 

16,8 

8.5 

8.7 

18.4 

16.4 

8.0 

MUGU (U - 48km) 

s 

16".4 

14.8 


10.0 

4.4 

11.5 

13 .’5 

20.0 ■ 

17.4 

• 10.5 

WSMR (U - 48KM) 

11.0 

19.5 

15.5 

15.1 

8.0 

8.8 

11.6 

13.4 

24.6 

18.8 

11.3 

BARKING SANDS (U - 48KM) 





14.3 



23.9 

26,5 

25.0 

23.5 

ANNUAL AMPLITUDES 












GREELY (U - 48KM) 





2l,3 

25.5 

32.9 

27,1 

i3.2 

■ 28.6 


WALLOPS (U - 48KM) 

49.8 


43.0 

A 

74.5 

60.0 

57.6 

61.5 

59.1 

65,8 

61.9 

MUGU (U - 48KM) 


53.9 

46.0 


63’. 6 

60.4 

49.1 

51.8 

46.8 

55.6 

53.4 

WHITE SANDS (U - 48KM) 

54.4 

50,3 

46,0 

58. ,3 

63.8 

56,7 

48.0 

49,0 

43,5 

48.4 

52,4 

BARKING SANDS (U 48KM) 




* 

42.4 



36,1 

■37.6 

•:27.9 

: 38.2 
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SEASONS winter 



SEASONS SPRING 



POWER <M?/SeC2) 

VAN 

N P.R.t . ’ 

POWER <w2/SbC2) 

VAR 

N P,*' 


KH « 


^0 

« 

00 

303 

534 

406 

236 

2B3 

377 

397 

305 

23/ 

1|H 

«97 

3 

159330 

« 

11 

26 

51 

58 

55 

42 

48 

71 

74 

90 

58 

61 

U 

1 351 



60 

301 

456 

33S 

243 

329 

401 

379 

272 

203 

104 

44J 

1 I 

135135 

« 

4 

41 

66 

69 

68 

55 

53 

74 

83 

96 

59 

94 

34 

9b 




59 

341 

411 

2R1 

237 

314 

J62 

336 

234 

1d7 

B1 

367 

30 

102102 


11 

70 

99 

97 

B6 

69 

66 

Bb 

92 

92 

SO 

117 

73 

51 

1 


o 

119 

434 

500 

J lr« 

20J 

20 7 

243 

26M 

206 

139 

71 

397 

61 

64 

64 

« 

24 

84 

116 

116 

92 

69 

73 

87 

84 

74 

30 

132 

U1 

21 

21 


•> 

16S 

4RS 

516 

319 

IBS 

149 

IbS 

191 

169 

119 

61 

35b 

102 

24 

24 

o 

34 

81 

192 

100 

82 

66 

6b 

71 

67 

60 

30 

10$ 

139 

21 

2J 


« 

IB4 

420 

402 

27b 

1R2 

122 

137 

155 

UR 

84 

40 

116 

115 

23 

21 

* - 

38 

70 

7b 

76 

75 

67 

59 

59 

57 

46 

21 

91 

153 

20 

20 


« 

I6P 

326 

287 

22S 

164 

92 

109 

126 

07 

52 

22 

223 

123 

23 

23 

« 

35 

55 

58 

67 

72 

b5 

57 

59 

54 

42 

20 

82 

164 

20 

20 


o 

145 

280 

259 

212 

152 

B1 

as 

94 

65 

43 

IR 

193 

130 

22 

2? 

4 

35 

4b 

50 

62 

68 

60 

55 

59 

52 

44 

23 

81 

172 

20 

20 


o 

130 

272 

26S 

227 

1 ?9 

n 

74 

81 

6? 

48 

24 

219 

113 

22 

22 

C 

4l 

54 

49 

55 

brt 

50 

48 

bl 

43 

39 

21 

7b 

182 

19 

19 


* 

116 

242 

273 

214 

los 

49 

62 

69 

55 

46 

24 

184 

138 

22 

22 

4 

49 

61 

49 

4B 

47 

36 

36 

37 

31 

30 

16 

62 

laa 

IS 

19 

40 


90 

197 

226 

176 

Hb 

46 

61 

b'i 

42 

31 

15 

l4Cl 

144 

22 

2? 

«• 

60 

59 

44 

42 

40 

30 

26 

2b 

25 

23 

U 

57 

189 

19 

19 



9? 

IR? 

19S 

140 

7b 

60 

66 

52 

37 

26 

9 

137 

145 

22 

22 


4B 

53 

37 

35 

33 

24 

22 

26 

2b 

20 

9 

48 

192 

19 

is 


» 

06 

176 

177 

119 

77 

70 

70 


39 

33 

12 

129 

14b 

21 

21 

• 

46 

51 

33 

26 

25 

17 

18 

25 

2b 

20 

B 

4b 

192 

19 

19 


O’ 

HS 

166 

164 

116 

B3 

72 

67 

45 

J4 

33 

14 

132 

143 

22 

22 

-» 

4S 

52 

34 

26 

22 

16 

17 

21 

23 

20 

B 

41 

1 94 

19 

19 


» 

RA 

16S 

153 

106 

79 

66 

52 

3b 

27 

2? 

12 

116 

14b 

21 

21 

« 

43 

51 

33 

25 

21 

16 

15 

16 

Id 

17 

8 

«»3 

201 

19 

19 

30 

» 

9? 

169 

147 

9S 

73 

S'y 

4J 

jy 

2b 

23 

10 

105 

146 

22 

22 

* 

42 

50 

32 

23 

20 

16 

13 

13 

16 

lb 

7 

3U 

204 

19 

19 


«o 

«e«« 


>40o< 

BaS040444 «««»«««« 












P£W* 

44 

2P 

14. 

T ll 

M.rt 

7,3 

6,3 

5.S 

4.9 

4.4 

4 

<UAYS> 



PhPlUO 

44 

22 

14,7 

11 

8,8 

7.3 

6.3 

5.5 

4.9 

4.4 

4 

(DAYS) 




StASDNs 

P0WE«» t«t/SFCrf3 


rCM e 


60 


6 

30 

3? 

31 

3j 

37 

SO 

65 

66 

54 



6 

2P 

34 

Jb 

is 

40 

44 

53 

54 

4J 


» 

5 

2? 

32 

35 

37 

35 

32 

34 

32 

2b 



3 

17 

2b 

26 

2b 

2*> 

24 

23 

21 

22 


o 

? 

IS 

23 

21 

iw 

21 

21 

14 

20 

22 

50 

« 

1 

12 

20 

17 

14 

16 

17 

17 

Ift 

Lb 


» 

2 

Q 

14 

12 

1 1 

11 

12 

U 

13 

13 



1 

7 

10 

10 

10 

11 

11 

l‘> 

10 

10 


« 

1 

S 

9 

9 

10 

11 

11 

10 

9 

8 


o 

1 

4 

7 

e 

d 

0 

9 

w 

7 

b 

40 

« 

1 

4 

5 

b 

b 

5 

5 

b 

4 

5 


* 

1 

4 

4 

4 

4 

4 

4 

3 

3 

4 


« 

1 

3 

5 

S 

5 

4 

4 

3 

3 

3 



1 

3 

4 

4 

4 

3 

3 

3 

3 

3 


o 

1 

2 

3 

J 

3 

3 

3 

3 

3 

3 

30 

• 

1 

? 

3 

3 

3 

2 

<r 

3 

3 

3 


PEP* 2? 14.7 11 7«3 6.3 b«b> 4«<) 4.4 


REASONS AUTUKIJ 



VAR 

N 

P.( 

<.e. 

. POWER 

(M2/5EC21 

2S 

bU 

23 

10bio» 

1 

10 

10b 

194 

177 

|9 

6U 

75 

49 

49 

« 

8 

sa 

157 

135 

12 

4i 

11b 

23 

23 

* 

1 

81 

106 

82 

12 

2/ 

134 

22 

2? 

4 

4 

67 

94 

85 

12 

29 

142 

22 

22 

» 

b 

bb 

107 

111 

9 

21 

145 

22 

22 

4 

17 

bl 

100 

100 

6 

16 

147 

22 

2? 

B 

28 

58 

76 

80 

S 

11 

151 

22 

22 

B 

f 

56 

6b 

69 

4 

14 

152 

22 

22 

« 

< 

50 

64 

69 

4 

V 

155 

22 

2? 

B 

2* 

47 

52 

52 

3 

6 

ISB 

21 

99 

c- 

31 

49 

43 

37 

2 

4 

158 

22 

22 

« 

33 

4R 

39 

32 

2 

6 

156 

22 

2? 

• 

29 

r 

32 

26 

2 

3 

154 

22 

22 

« 

27 

3b 

2/ 

23 

2 

4 

154 

22 

22 

o 

25 

33 

24 

21 

2 

3 

154 

22 
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Table A-1, Analysis of the high frequency variability of the wind at 
Greely; (a) Zonal, (b) Meridional. 
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22 
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0 
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89 
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74 
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87 
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Table A-2. Analysis "of the 'high frequency variability of t'he 
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Table A-3. Analysis of the high frequency variability of the 
wind at Wallops, (a) Zonal, (b) Meridional. 
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30 196 

0 

2 

14 

20 

21 

20 

- 19 

21 

21 

20 

JB 

9 

26 232 

0 

1 

11 

19 

20 

18 

16 

17 

17 

17 

16 

a 

21 253 

0 

1 

10 

17 

16 

15 

14 

14 

14 

14 

IS 

6 

20 263 

0 

1 

ID 

15 

14 

12 

11 

11 

12 

13 

12 

6 

16 267 

0 

1 

8 

12 

10 

8 

9 

10 

11 

11 

10 

s 

13 268 

0 

1 

6 

9 

7 

7 

9 

10 

10 

9 

8 

4 

11 268 

0 

0 

4 

6 

6 

6 

8 

9 

9 

• 6 

6 

3 

8 269 


9 266 
6 266 
4 265 
4 263 
J 261 


PER. 44 22 U.7 11 a. 8 7*3 6.3 5*5 4*9 4,4 


109109 
70160 
19 19 
18 18 
J7 17 
17 17 
16 16 
16 16 
16 16 
16 16 
16 16 
16 16 
16 16 
16 16 
16 16 
16 16 


VAR P.R.E* 


24 66 

23 67 


17 

11 

9 

9 


8 33 


30 
29 
25 
20 
17 . 
15 

n 

6 

4 


4 (OAYSJ 


83 

75 

S7 

59* 

67 

63 

75 

75 

34 

110 ' 15 

120] 

78 

69 

> 54 

49 

5U 

48 

63 

•69 

32 

02 45 

80 

v69 

60 

.50 

37 

28 

30 

50 

58 

27 

♦ 66 99 

22 

63 

S3 

44 

* 34 

24 

28 

43 

47 

22 

59 146 


$4 

45 

39 

32 

25 

28 

35 

37 

10 

50 170 

19 

44 

36 

36 

31 

24 

27 

31 

•27 

13 

45* 185 

18 

36 

31 

34 

'26 

21 

25 

25 

18 

. 7 

36 199 

18 

3% 

29 

.29 

23 

16 

'20 

20 

14 

6 

31 205 


34 

27 

22 

17 

12 

14 

17 

15 

7 

27 207 

10 

32 

26 

18 

13 

10 

12 

16 

IS 

7 

25 208 

18 

28 

22 

IS 

12 

10 

12 

15 

13 

5 

22 208 

'18 

22 

17 

13 

’ 12 

10 

12 

12 

10 

5 

18 208 

18 

17 

13 

11 

1 ll 

9 

9 

9 

9 

5 

IS 207 

18 

13 • 

11 

• 10 

9 

7 

6 

7' 

7 

4 

207 

16 

9 

10 

9 

8 

7 

6 

S 

6 

3 

9 208 

18 

7 

9 

8 

8 

7 

6 

5 - 

S 

3 

/ 207 

18 

14.7 .11 

8^8 

7.3 6 

.3 

0B0* 

5.5 

4.9 4 

• ••4 

• 4 

»00 

4 

(OAYSl 



Table A-4. . Analysis of the high frequency variability- of the 
wind at Pt, Mugu. ' (a) Zonal, (b) Meridional. 



57 


^£aSOMs WINTfH 
POWEft <H?/SEC^} 


KH o 


60 


102 

369 

464 

34d 

251 

209 

302 

171 

99 

115 


« 

92 

36P 

479 

366 

263 

29e» 

300 

103 

113 

117 


« 

97 

372 

402 

303 

273 

203 

205 

170 

132 

131 



UP 

30A 

466 

3?a 

267 

262 

254 

152 

120 

132 


9 

127 

399 

465 

363 

252 

246 

233 

139 

99 

113 

50 

« 

lAl 

A23 

477 

3*5 

239 

217 

192 

12B 

101 

107 


ft 

156 

A31 

476 

34* 

235 

195 

15B 

110 

110 

111 


O 

163 

4Q4 

426 

320 

227 

176 

142 

115 

105 

107 


ft 

170 

360 

36? 

263 

201 

159 

131 

103 

04 

02 


ft 

173 

369 

356 

240 

164 

141 

117 

07 

67 

55 

AO 

ft- 

159 

339 

326 

205 

134 

lln 

99 

01 

65 

41) 


ft 

139 

293 

290 

105 

117 

94 

74 

60 

62 

36 


ft 

120 

2AS 

250 

166 

100 

71 

51 

49 

50 

33 


ft 

Qft 

19? 

199 

13S 

7h 

52 

30 

3* 

34 

27 


O' 

76 

1A2 

US 

99 

56 

35 

20 

27 

25 

20 

30 

» 

65 

117 

119 

03 

4b 

26 

22 

25 

23 

17 


PE»), A4 2^ 14,7 11 a«S 7«J 6*3 b«b 4.Q 4«4 


b£4S0H= SPklNC 



VAK 

N 

P.fl 

I.E. 

HOXER 

(H2/SEC2) 


60 

322 

90 

31 

31 

• 

ft 

26 

98 

117 

99 

57 

407 

126 

20 

20 

♦ 

29 

08 

106 

92 

57 

359 

154 

19 

19 

9 

33 

74 

06 

70 

60 

305 

171 

19 

19 

ft 

35 

68 

73 

70 

59 

33b 

102 

10 

10 

ft 

34 

64 

68 

69 

56 

349 

190 

10 

10 

ft 

31 

61 

60 

67 

59 

34U 

195 

10 

in 

ft 

31 

63 

O0 

61 

61 

329 

190 

10 

10 

ft 

30 

64 

70 

63 

40 

203 

202 

10 

10 

ft 

31 

65 

74 

69 

29 

263 

205 

10 

10 

« 

33 

63 

68 

59 

16 

22* 

209 

10 

10 

ft 

33 

59 

50 

45 

a 

199 

213 

18 

10 

• 

30 

57 

55 

40 

a 

167 

215 

18 

16 


26 

S2 

52 

39 

10 

130 

213 

18 

IK 

9 

20 

41 

*2 

36 

6 

9* 

212 

18 

10 

9 

16 

32 

34 

31 

6 

65 

214 

18 

IH 

• 

15 

32 

33 

29 

*99 

4 

(DAYS) 



PEHIUD 

aftftf 

44 

^ftfte( 

22 

^e«e•{ 

14.7 

aftfti 

a 


VAR N P.P«E. 


00 

56 

57 

97 

126 

110 

51 

150 

140 

19 

92 

69 

50 

60 

72 

98 

103 

54 

109 

100 

10 

18 

57 

46 

45 

42 

50 

06 

53 

91 

209 

10 

18 

6l 

S3 

46 

35 

39 

63 

42 

83 

233 

1 7 

17 

67 

56 

41 

29 

31 

49 

33 

77 

251 

17 

17 

65 

52 

36 

27 

30 

49 

33 

73 

260 

17 

17 

54 

46 

39 

31 

32 

47 

29 

73 

268 

17 

17 

53 

47 

44 

36 

34 

36 

10 

67 

271 

17 

17 

59 

49 

41 

31 

20 

26 

11 

76 

273 

17 

17 

51 

42 

32 

24 

22 

20 

10 

50 

274 

16 

16 

38 

32 

26 

21 

20 

10 

9 

52 

276 

16 

16 

33 

29 

24 

19 

19 

17 

6 

46 

201 

16 

16 

31 

27 

22 

10 

18 

16 

7 

46 

284 

16 

16 

30 

24 

19 

17 

17 

15 

7 

30 

286 

16 

16 

26 

19 

15 

14 

13 

12 

6 

31 

287 

16 

16 

23 

16 

13 

12 

a 

10 

5 

26 

287 

16 

16 

8.0 

>••«« 

7.3 

»ftftfti 

6.3 

iftftft. 

5.5 

tftftO. 

4.9 

4.4 

4 

{DAY5) 




SeASON= SUMMER 
POUER ^«^/SEC^J 


VAK N P,R*£. 


SEASOM^ AUTUMM 
HOi*ER CH2/SEC2) 


VAK « P.R»E* 


60 

ft 

? 

6S 

112 122 

1U4 

95 

93 

9* 

107 

65 

31 

150 

137 

19 

19 

ft 

14 

114 

155 

157 

ai 

149 

lU 

116 

97 

70 

45 

160 

127 

19 

19 


ft 

6 

66 

102 102 

92 

69 

07 

91 

99 

75 

20 

105 

1S8 

17 

17 

ft 

15 

ll)6 

149 

153 

160 

130 

107 

112 

loo 

62 

45 

170 

180 

18 

10 


ft 

9 

57 

02 7* 

70 

79 

00 

65 

05 

68 

31 

110 

229 

17 

a 

ft 

20 

104 

152 

156 

146 

117 

94 

107 

106 

90 

45 

146 

213 

16 

18 


ft 

7 

3P 

59 56 

5* 

66 

74 

75 

69 

64 

35 

01 

253 

17 

a 

ft 

30 

117 

168 

164 

140 

109 

91 

107 

112 

96 

47 

181 

227 

17 

17 


ft 

S 

27 

40 *2 

42 

ss 

66 

64 

53 

50 

20 

65 

268 

16 

16 

ft 

35 

120 

172 

152 

114 

90 

93 

97 

93 

as 

45 

162 

236 

17 

17 

SO 

ft 

3 

2? 

34 3H 

39 

45 

54 

52 

41 

38 

22 

55 

279 

16 

16 

ft 

38 

110 

15U 

119 

75 

77 

0* 

75 

61 

56 

32 

U3 

252 

17 

17 



2 

19 

32 37 

37 

36 

*0 

39 

32 

29 

16 

43 

206 

16 

16 

e 

43 

9b 

112 

06 

50 

64 

67 

57 

45 

39 

22 

92 

264 

17 

17 



3 

IP 

29 33 

33 

30 

31 

32 

29 

24 

12 

37 

260 

16 

16 

ft 

*»* 

89 

9J 

72 

53 

55 

55 

48 

42 

37 

19 

87 

267 

16 

16 



2 

16 

25 20 

28 

25 

24 

29 

29 

24 

a 

35 

290 

16 

16 


ft 

07 

85 

69 

54 

51 

49 

43 

37 

31 

16 

61 

271 

16 

16 



3 

17 

25 28 

29 

24 

19 

23 

26 

22 

10 

28 

294 

16 

16 

ft 

4 ' 

82 

7b 

69 

50 

49 

44 

39 

31 

24 

12 

75 

275 

16 

16 

40 


3 

17 

24 28 

31 

26 

19 

19 

20 

18 

9 

35 

300 

16 

16 

ft 

46 

04 

79 

60 

52 

37 

34 

33 

29 

25 

13 

72 

276 

16 

16 



2 

12 

10 22 

24 

20 

17 

15 

15 

14 

7 

21 

303 

16 

16 

• 

48 

92 

04 

63 

*0 

25 

22 

25 

26 

27 

14 

69 

277 

16 

16 



1 

7 

13 a 

15 

13 

14 

13 

11 

11 

6 

15 

304 

16 

16 


44 


76 

54 

33 

20 

16 

20 

26 

24 

a 

59 

27? 

16 

16 



0 

6 

13 16 

13 

a 

13 

12 

11 

10 

S 

17 

303 

16 

16 

ft 

36 

65 

50 

42 

29 

20 

17 

19 

21 

19 

9 

47 

277 

16 

16 



1 

6 

12 14 

12 

10 

10 

9 

9 

8 

4 

12 

302 

16 

16 

• 

29 

51 

43 

33 

25 

18 

17 

17 

a 

17 

10 

39 

277 

16 

16 

30 

ft 

1 

6 

11 12 

10 

8 

7 

/ 

7 

7 

3 

10 

304 

16 

16 

ft 

ft 

27 

46 

30 

29 

22 

16 

15 

15 

15 

17 

10 

33 

270 

16 

16 

ft. 

PER- 

44 

«ft«< 

22 

14.7 a 

• ftfti 

0.0 

ftft«< 

7.3 

**« 

6.3 

5,5 

.ftftfti 

4.9 

Oftftft 

4.* 

»»« 

4 

(DAYS) 



• ft 

PERIOD 

tftft 

44 

.••ft 

22 

• ftO 

U. 

7 a 

«ft*< 

0.8 

7.3 

6.3 

5-5 

4.9 

•••ft 

4.4 

• •• 

4 

(DAYS) 




>) 


so 


AO 


30 


SEASONS VINTEH 
POUEP (H2/SEC«!J 



1 

107 

150 

162 

156 

132 

as 

117 

115 

114 


a 

118 

170 

167 

152 

122 

90 

101 

107 

99 


29 

131 

167 

151 

130 

100 

71 

01 

100 

90 


43 

138 

140 

125 

107 

79 

6l 

76 

95 

98 


44 

137 

146 

129 

108 

74 

65 

74 

62 

65 


42 

133 

151 

138 

110 

77 

71 

76 

77 

68 


44 

129 

140 

135 

106 

77 

71 

73 

74 

59 


37 

109 

130 

122 

96 

73 

60 

64 

60 

48 


2P 

90 

m 

99 

77 

61 

60 

56 

48 

41 


25 

78 

90 

74 

61 

50 

49 

48 

42 

39 


21 

58 

64 

56 

S3 

43 

35 

33 

33 

34 


14 

40 

50 

51 

40 

36 

25 

22 

26 

20 

ft 

0 

30 

4A 

47 

40 

32 

23 

21 

23 

23 


5 

22 

33 

35 

30 

26 

22 

20 

19 

20 


5 

16 

23 

23 

21 

20 

a 

13 

12 

15 


6 

15 

19 

18 

16 

17 

15 

10 

0 

12 


S£A50N= SPRING 



VAR 

N 

P.fi 

• E. 

POWER 

(M2/5EC21 

61 

168 

90 

31 

31 

« 

ft 

2 

31 

56 

66 

49 

176 

126 

20 

20 

ft 

1 

29 

51 

57 

44 

146 

150 

19 

19 

ft 

0 

26 

42 

43 

53 

U6 

165 

19 

19 

ft 

2 

23 

36 

36 

47 

138 

177 

19 

19 

• 

6 

22 

29 

30 

32 

136 

185 

18 

18 

ft 

6 

21 

27 

29 

24 

138 

189 

18 

18 

ft 

6 

22 

30 

33 

20 

la 

192 

10 

10 

ft 

6 

22 

32 

36 

21 

94 

193 

IS 

18 

ft 

S 

21 

32 

34 

21 

84 

195 

18 

le 

ft 

.5 

20 

29 

28 

19 

62 

200 

18 

10 

ft 

6 

19 

25 

23 

14 

46 

205 

16 

10 

ft 

6 

16 

19 

20 

12 

45 

209 

18 

10 

ft 

6 

13 

15 

16 

a 

34 

206 

IB 

16 

ft 

5 

12 

13 

13 

10 

23 

205 

18 

IB 

ft 

3 

10 

12 

a 

0 

17 

208 

18 

10 

ft 

2 

6 

a 

10 


VAR W P*R.E. 


63 

70 

77 

70 

75 

63 

30 

95 

140 

19 

19 

55 

63 

70 

71 

70 

66 

34 

76 

180 

18 

18 

44 

51 

57 

57 

59 

62 

33 

68 

209 

18 

18 

2B 

44 

AS 

45 

46 

44 

23 

52 

233 

17 

17 

33 

40 

40 

37 

37 

38 

21 

41 

251 

17 

17 

33 

40 

38 

34 

36 

42 

25 

51 

260. 

a 

17 

37 

40 

35 

29 

32 

39 

23 

44 

268 

17 

17 

36 

36 

32 

27 

29 

33 

19 

45 

271 

17 

17 

32 

30 

29 

27 

27 

26 

14 

39 

273 

17 

17 

26 

27 

26 

24 

20 

a 

8 

32 

274 

16 

16 

22 

23 

23 

19 

IS 

12 

6 

26 

276 

16 

16 

20 

20 

19 

IS 

12 

12 

6 

23 

281 

16 

16 

16 

16 

16 

13 

12 

13 

7 

19 

284 

16 

16 

12 

12 

14 

13 

12 

12 

7 

18 

286 

16 

16 

9 

10 

12 

10 

9 

ID 

6 

13 

287 

16 

16 

9 

10 

10 

8 

6 

8 

5 

11 

207 

16 

16 


PER. 4A 2Z 1A.7 II 8.8 7.3 6,3 5.5 A. 9 A.A A (OAYSI 


PERIOD AA 22 1A.7 U 0.8 7.3 6.3 5.S A. 9 A.A A <OAYS) 


SEASONS SUMMER 
POWER {M2/SEC2> 


KH • 


60 * 

1 

41 

65 

76 

76 

69 

61 

52 

61 

67 

0 

1 

39 

65 

72 

66 

58 

SS 

50 

56 

62 

• 

0 

33 

59 

bS 

55 

48 

51 

4b 

46 

50 

ft 

0 

26 

47 

54 

53 

52 

SO 

44 

39 

38 

ft 

I 

22 

37 

J9 

43 

47 

42 

38 

36 

32 

50 « 

2 

21 

31 

29 

30 

36 

35 

33 

33 

28 

ft 

3 

19 

30 

27 

25 

29 

30 

28 

27 

24 

ft 

2 

16 

26 

26 

24 

23 

22 

22 

22 

20 

ft 

2 

12 

18 

22 

22 

10 

a 

19 

19 

a 

ft 

2 

0 

13 

16 

17 

15 

14 

16 

10 

16 

40 • 

1 

7 

11 

13 

14 

13 

12 

12 

14 

14 

ft 

1 

7 

10 

a 

12 

12 

a 

9 

10 

a 

ft 

1 

6 

9 

9 

9 

10 

9 

9 

9 

9 

ft 

0 

4 

7 

6 

8 

8 

0 

8 

8 

8 

ft 

0 

3 

6 

7 

7 

7 

7 

7 

7 

6 

30 * 

ft 

0 

3 

6 

7 

7 

7 

7 

7 

6 

6 

ftft 

PER. 

ftftftft 

44 

• ••4 

22 

^ftftftft< 

U.7 

»•»< 

a 

»•••! 

0.8 

Kftft«< 

7.3 

>««•« 

6«3 

•ftfti 

5.5 

^ftftfti 

4*9 

(ftftfti 

4.4 


SEASONS AUTUMN 



VAR 

N 

P.R 


POWER 

(M2/SEC2J 


33 

88 

137 

19 

19 

ft 

ft 

3 

46 

79 

92 

31 

79 

lae 

17 

17 

ft 

1 

44 

71 

68 

26 

63 

229 

a 

17 

ft 

0 

45 

62 

77 

20 

64 

253 

17 

17 

ft 

a 

51 

65 

63 

16 

46 

266 

16 

16 

* 

12 

53 

65 

56 

13 

41 

279 

16 

16 

ft 

9 

49 

64 

5S 

a 

36 

286 

16 

16 

ft 

7 

44 

61 

57 

10 

30 

288 

16 

16 

ft 

6 

40 

52 

48 

8 

24 

290 

16 

16 

ft 

6 

34 

44 

42 

7 

19 

294 

16 

16 

ft 

4 

32 

41 

37 

7 

16 

299 

16 

16 

ft 

5 

29 

36 

31 

5 

13 

303 

16 

16 

ft 

5 

23 

26 

23 

5 

12 

304 

16 

16 

ft 

5 

18 

22 

19 

4 

9 

303 

16 

16 

ft 

3 

15 

20 

18 

3 

8 

302 

16 

16 

ft 

2 

U 

15 

16 

3 

0 

304 

16 

16 

ft 

1 

9 

13 

14 

^«• 

4 

(DAYS) 



• «< 

PERIOD 

iftfti 

44 

^es•l 

22 

14.7 

a 









VAR 

N 

P.R 

:*£. 

94 

84 

78 

61 

72 

77 

49 

as 

127 

19 

19 

91 

80 

73 

73 

69 

74 

44 

94 

179 

IB 

16 

83 

74 

63 

57 

59 

61 

33 

92 

213 

16 

10 

66 

65 

52 

44 

44 

41 

21 

70 

226 

17 

17 

b5 

54 

43 

36 

36 

30 

13 

60 

234 

17 

17 

b6 

50 

38 

33 

35 

29 

a 

62 

249 

17 

17 

53 

49 

37 

30 

31 

26 

11 

57 

261 

17 

17 

47 

44 

34 

29 

29 

25 

a 

50 

264 

16 

16 

43 

38 

26 
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Table A-5. Analysis of the high frequency variability of the 
wind at XJhite Sands. fa) Zonal, (b) Meridional. 
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Table A-6. Analysis of the high frequency variability -of the 
wind at Kennedy. (a) Zonal, .(b) IJeridional. 
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Table A-7. Analysis of the high frequency variability of the 
wind at Hawaii, (a) Zonal, (b) Meridional. 
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Table A- 8. 


Analysis 6f the high frequency variability* "bf the 
wind at Ascension. 



